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APPARATUS AND METHOD FOR DETERMINING THE 
ACTIVE DOPANT PROFILE IN A SEMICONDUCTOR WAFER 

5 DISCUSSION OF THE RELATED ART 

In the processing of a semiconductor wafer to form integrated circuits, 
charged atoms or molecules are directly introduced into the wafer in a process 
called ion implantation. Ion implantation normally causes damage to the lattice 
structure of the wafer, and to remove the damage, the wafer is normally annealed at 

10 an elevated temperature, typically 600°C to 1 lOO^C. This anneal also causes 

implanted atoms to move from interstitial sites to substitutional sites in the crystal 
lattice (an atom must be in a substitutional site to be electrically active). Prior to 
annealing, material properties at the surface of the wafer may be measured, 
specifically by using the damage caused by ion implantation. 

15 For example, U.S. Patent 4,579,463 granted to Rosencwaig et al. (that is 

incorporated herein by reference in its entirety) describes a method for measuring a 
change in reflectance caused by a periodic change in temperature of a wafer's 
surface (see column 1, lines 7-16). Specifically, the method uses "thermal waves 
[that] are created by generating a periodic localized heating at a spot on the surface 

20 of a sample'' (column 3, lines 54-56) with "a radiation probe beam ... directed on a 
portion of the periodically heated area on the sample surface," and the method 
"measur[es] the intensity variations of the reflected radiation probe beam resulting 
from the periodic heating'' (column 3, lines 52-66), 

As another example, U.S. Patent 4,854,710 to Opsal et al. (also 

2 5 incorporated herein by reference in its entirety) describes a method wherein "the 

density variations of a diffusing electron-hole plasma are monitored to yield 
information about features in a semiconductor" (column 1, lines 61-63). 
Specifically, Opsal et al. state that "changes in the index of reflection, due to the 
variations in plasma density, can be detected by reflecting a probe beam off the 

3 0 surface of the sample within the area which has been excited" (column 2, lines 23- 

3 1) as described in "Picosecond EUipsometry of Transient Electron-Hole Plasmas 
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in Germanium," by D. H. Auston et al.. Physical Review Letters, Vol. 32, No. 20, 
May 20, 1974. Opsal et al. further state (in column 5, lines 25-3 1 of U.S. Patent 
4,854,710): "The radiation probe will undergo changes in both intensity and phase. 
In the preferred embodiment, the changes in intensity, caused by changes in 
5 reflectivity of the sample, are monitored using a photodetector. It is possible to 
detect changes in phase through interferometric techniques or by monitoring the 
periodic angular deflections of the probe beam." 

A brochure entitled "TP-500: The next generation ion implant monitor" 
dated April, 1996 published by Therma-Wave, Inc., 1250 Reliance Way, Fremont, 

10 CA 94539, describes a measurement device TP-500 that requires "no post-implant 
processing" (column 1, lines 6-7, page 2) and that "measures lattice damage" 
(colunm 2, line 32, page 2). The TP-500 includes "[t]wo low-power lasers [that] 
provide a modulated reflectance signal that measures the subsurface damage to the 
silicon lattice created by implantation. As the dose increases, so does the damage 

15 and the strength of the TW signal. This non-contact technique has no harnifiil 
effect on production wafers" (columns 1 and 2 on page 2). According to the 
brochure, TP-500 can also be used after annealing, specifically to "optimize ... 
system for annealing uniformity and assure good repeatability" (see bottom of 
colimm 2, on page 4). 

20 

SUMMARY 

An apparaUis and method in accordance with the invention stimulate a • 
region of a semiconductor wafer (also called "semiconductor substrate") that 
originally has a first number of charge carriers, so that there are a second number of 

2 5 charge carriers during the stimulation. The stimulation can be accomplished in any 

number of ways, including e.g. by use of a beam of electromagnetic radiation or by 
a beam of electrons. The apparatus and method use a measurement device (such as 
an interferometer in one embodiment) to obtain a measured value of a signal that is 
affected by the stimulation. In one embodiment, the affected signal is a probe 

3 0 beam that is reflected by the charge carriers, although other signals can be used in 

other embodiments. 
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The apparatus and method also operate a simulator (e.g. a personal 
computer programmed with simulation software) to generate a simulated value for 
the measured signal. The simulated value is based on: (i) conditions present during 
stimulation (as described above) and (ii) a predetermined profile of the 
5 concentration of active dopants in the region under stimulation. If the measured 
value matches the simulated value, then the predetermined profile used in 
simulation is used as a measure of the profile of active dopants in the region. The 
simulation may be repeated with a number of such predetermined profiles. 
In one implementation, the simulations are repeated (prior to the 

10 stimulation) to obtain a set of such profiles, and the corresponding simulated 
values are used later to obtain a measure of the profile of active dopants in the 
region, e.g. by finding the closest simulated value to the measured value. In 
another implementation, one or more simulations are repeated after the stimulation 
only in case there is no match, until the simulated value and the measured value 

15 differ by less than a predetermined amount (e.g. less than 1%), and the 

corresponding predetermined profile is used as a measure of the profile of active 
dopants in the region. 

The measured profile of active dopants can be used in a number of ways. 
In one embodiment, the measured profile is used to determine junction depth that is 

2 0 compared with specifications for acceptability of the wafer. If the junction depth 
falls within the specifications, the wafer is processed fiirther (e.g. in a wafer 
processing unit to form another layer on the substrate, or in an annealer for heat 
treatment of the substrate), and otherwise the substrate is identified as unacceptable 
and placed in a bin of rejected substrates. 

25 In one embodiment, the apparatus and method creates charge carriers in a 

region of the semiconductor material (also called "carrier creation region") in a 
concentration that changes in a periodic manner (also called "modulation") only 
wdth respect to time. Thereafter, the apparatus and method determine the number 
of charge carriers created in the carrier creation region by (I) measuring an 

30 interference signal obtained by interference between a reference beam and a portion 
of a probe beam that is reflected by the charge carriers, and (2) comparing the 



-3- 



wo 00/57159 ^ ^ PCT/USOO/07357 



measurement with predetemiined data (e.g. in a graph of such measurements 
plotted against junction depth). 

Charge carriers that are created as described above (also called "excess 
carriers") are in excess of a number of charge carriers (also called "background 
5 carriers") that are normally present in the semiconductor material m the absence of 
illumination. The concentration of excess carriers is modulated in time at a 
frequency that is maintained sufficiently small to ensure that the variation in 
concentration is aperiodic (i.e. not oscillatory, e.g. decays exponentially or 
according to a monotonic flinction). Specifically, a profile of excess carrier 

1 0 concentration that is devoid of a wave (along radial distance) is created as 

described herein when a:t least a majority (i.e. greater than 50%) of the charge 
carriers that move out of the carrier creation region do so due to diffusion. 

Such a temporal modulation under diffusive conditions (also called 
"diffusive modulation") is used to measure an interference signal, (for example, the 

15 phase and amplitude are both measured). The measurement is used to determine 
(e.g. by looking up a graph or a table) one or more properties (also called 
"semiconductor properties") of the semiconductor material (such as junction 
depth). The concentration of excess carriers as a function of depth from the front 
surface of the semiconductor material, when measured as described herein, can 

2 0 also be used to determine the concentration of active dopants in the semiconductor 
material. Specifically, a profile of excess carrier concentration in depth is a 
function of the depth profile of the electric field that results from the active dopants 
(that form the doped semiconductor material). 

An increase in excess carriers as a function of depth causes a corresponding 

2 5 increase in an index of refraction of the semiconductor material. Therefore, a laser 

beam (called "probe beam") shone on the semiconductor material is reflected back 
(by both background carriers and by excess carriers, but only the reflection by the 
excess carriers varies periodically at the modulation frequency), and a signal at the 
modulation frequency generated by interference between the reflected portion and a 

3 0 reference beam is measured as described herein. Various properties of the 

interference signal (such as amplitude and phase) are functions of the depth at 
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which the reflection occurs, and can be measured to determine the depth of the 
junction. Note that as used herein, a junction is the boundary of any doped region 
(irrespective of whether the doping is a p-type dopant into an n-type substrate, a p- 
type dopant into a p-type substrate, or vice- versa). 
5 A first embodiment (also called "front surface embodiment") measures the 

intensity of an interference signal that is obtained by interfering the reflected 
portion of the probe beam with a reference beam formed by another portion of the 
probe beam (this portion hereinafl:er being called "front surface beam") that is 
reflected by the front surface. In one variant of the front surface embodiment, a 

10 laser is used to generate another beam (called "generation beam") that is used to 
generate the excess carriers. The generation beam's intensity is modulated at a 
fixed frequency that is sufficiently low to ensure that the phase of the variation of 
the concentration of excess carriers is the same as (e.g, to vsdthin 10%) the phase of 
the generation beam over a diffusion length (wherein diffiision length is the length 

15 over which the excess charge carrier concentration decays to 1/e). Therefore, the 
excess carrier concentration changes approximately synchronously with the change 
in intensity of the generation beam. This condition ensures that the excess carrier 
distribution is primarily due to diffusion that can be modeled by a non-wave 
solution (rather than by a wave solution). 

20 In this variant, an interferometer measures the amplitude and phase of such 

an interference signal, as a function of the generation beam's power and 
modulation, and these measurements are used to determine the concentration of 
excess carriers. Variation in time of the excess carrier concentration as described 
above edlows the interferometer to use a lock-in amplifier to measure the reflected 

2 5 portion of the probe beam with an acctiracy not possible when the excess carrier 

concentration is fixed. 

In one implementation, a number of graphs relating the interference signal 
measurement to the junction depth and to the power of generation beam are 
determined (either by simulation or empirically). Thereafter, for a given wafer 

3 0 (also called "production wafer"), measurements (also called "interference 

measurements") of the interference signal for different powers of the generation 
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beam are performed, and the results are compared to one or more of the just- 
described graphs, thereby to determine a graph that indicates the junction depth. 
Specifically, predetermined graphs are generated in the following manner for a 
number of dopant profiles that approximate an expected dopant profile of the 
production wafer. 

Any method or device can be used to generate dopant profiles that are 
provided as input to a simulator (that may be a programmed computer executing a 
simulation program) for generation of the predetermined graphs. For example, 
spreading resistance profiles can be obtained on wafers (also called "reference 
wafers") that have been processed under known conditions and have known 
properties. Alternatively, dopant profiles can be simulated using commercially 
available simulators (that assume movement of charge carriers from the carrier 
creation region by diffusion). 

Next, for a given dopant profile, a profile of the excess carrier concentration 
as a ftmction of depth is determined using a simulator, for each of a number of 
powers of the generation beam. Next, a derivative of excess carrier profile as a 
function of depth z &om the fi-ont surface is multiplied by cos(2kn2), wherein 
k=2 ;r / A , with A being the wavelength of the probe beam, and n being the index 
of refi-action of silicon. The product of multiplication is integrated with respect to 
depth and multiplied by one or more constant factors (that are related to known 
physical constants and to calibration of the measurement system) to determine a 
simulated value of the interference measurement. The simulated value of the 
interference measurement is thereafter plotted on a graph as a function of depth z 
for a selected generation beam power. The just-described acts are repeated to 
obtain graphs for other fixed values (e.g. two additional values) of the generation 
beam power. Additional such graphs are generated for different dopant profiles. 

After such graphs are available, interference measurements on a production 
wafer at the selected generation beam power are used to look up the graphs to 
determine junction depth. The look up can be repeated for different interference 
measurements obtained by using different powers of the generation beam, to 
eliminate ambiguity that may result fi-om two wafers having different junction 
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depths but same measurements (as may occur, e.g. when changes in the two 
numbers being multiplied, namely (1) the derivative and (2) the cosine function (as 
described above) compensate for each other in the two wafers). A predetermined 
dopant profile having the same junction depth as that obtained by look up is 
5 thereafter used as the profile of active dopants present in the production wafer. 

Measurement of the phase and amplitude of an interference signal as 
described herein is a significant aspect of one implementation. One or more such 
measurements provide a measure of a property of the semiconductor material (or a 
process condition) during wafer fabrication. In another implementation (also 

10 called "scanning implementation"), a number of such measurements are performed 
at different locations on a wafer (while the generation beam's power is maintained 
constant). Any change in such measurements indicates a correspondihg change in 
the concentration of active dopants (at a predetermined depth from the front 
surface). Therefore, such interference measurements (from which active dopant 

15 profile is determined) are preferably (but not necessarily) monitored in one variant 
of the invention during wafer fabrication, to control a process step (e.g. to control 
annealing temperature of a wafer that has been ion implanted) used in fabricating 
the wafer. 

When the junction depth and junction profile are measured directly on the 
2 0 wafer undergoing fabrication (also called "patterned wafer" or "annealed wafer" 

depending on the stage of fabrication), a measurement as described herein 

increases yield, as compared to an off-line measurement of a test wafer's properties. 

Moreover, such a measurement avoids the prior art cost of the test wafer itself. 

Such measurements are performed in one embodiment after annealing a production 
25 wafer to activate the dopants, thereby to obtain a measvire that is more indicative of 

the electrical behavior of the devices being fabricated than a property that is 

measured prior to annealing (as described in U.S. Patent 4,854,710). 

In a second embodiment (also called "phase embodiment"), instead of the 

above-described interference signal, another interference signal is generated by 
30 interference between the reflected portion of the probe beam (described above) and 

another reference beam (hereinafter "variable phase beam") having a phase that can 
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be changed independent of the phase of the probe beam. A phase difference 
(detected using, e.g. a phase detector) between two interference signals indicates 
the junction depth, wherein a first interference signal is obtained by interference of 

(1) the variable phase beam and (2) the front surface beam that is described above 
5 as the portion of probe beam reflected by the front surface, and a second 

interference signal is obtained by interference of (1) the variable phase beam and 

(2) the reflected portion of the probe beam. 

In the second embodiment, the probe beam is coherent (i.e. of single 
chrominance, e.g. single wavelength) in addition to being polarized, so that 
10 interference with the variable phase beam can happen. Use of a reference beam in 
the second embodiment that is independent of the probe beam provides an increase 
in sensitivity of the measurement of material properties over the first embodiment, 
because of increased sensitivity of a phase detector used in the second embodiment 
to measure the interference signal. Use of an independent reference beam also 
allows absolute measurement of the junction depth as a fraction of the wavelength 
in the semiconductor. 



15 



20 



BRIEF DESC RIPTION OF THE PR A WTNOR 

FIG. 1 A illustrates, in a high level block diagram, a system including an 
apparatus (called "active dopant profiler") in accordance with the invention. 

FIGs. IB and IC illustrate, in high level flow charts, a method performed 
by the apparatus of FIG. 1 A. 

FIG. ID illustrates, in a graph, the temporal modulation of charge carriers 
by one embodiment of the active dopant profiler of FIG. lA, without creation of a 
2 5 wave of charge carriers. 

FIG. IE illustrates, in a cross-sectional view of the semiconductor, use of a 
probe beam, a generation beam, and an optional reference beam used by the active 
dopant profiler of FIG. 1 A in various embodiments described herein. 

FIG. IF illustrates, in another cross-sectional view, modeling of excess 
30 carriers in layers in the semiconductor. 
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FIG. IG illustrates, in a graph, the concentration of excess charge carriers 
as a function of depth Z from front surface 153, 

FIG. 2A illustrates, in a flowchart, the acts performed by the profiler of 
FIG. 1 A in one implementation. 

FIG. 2B illustrates, in a flow chart, creation of charge carriers in act 230 
(illustrated in FIG. 2 A) performed by the profiler of FIG. 1 A. 

FIG. 2C illustrates, in a flow chart, use of the measurements in optional act 
250 (illustrated in FIG. 2 A) performed by the profiler of FIG. 1 A. 

FIGs. 2D-2F illustrate various graphs containing the predetermined data 
that are used in one embodiment of act 210 illustrated in FIG, 2A. 

FIG. 3 illustrates, in a flowchart, generation of predetermined data in an 
optional act 210 (illustrated in FIG. 2 A) performed by the profiler of FIG. 1 A. 

FIG. 4A illustrates, in a graph, concentration of active dopants as a function 
of depth from the front surface of the semiconductor material (also called "active 
doping profile") for use in generation of predetermined data, obtained by 
performance of act 2 1 1 illustrated in FIG. 2B. 

FIG. 4B illustrates, in a graph, a number of profiles of excess carriers (i.e. 
the concentration of excess carriers as a function of depth) for the dopant profile of 
FIG. 4 A, obtained by performance of act 241 illustrated in FIG. 2 A. 

FIG. 4C illustrates, in a graph, a simulated value of an interference 
measurement as a function of the intensity (also called "power*') of the generation 
beam, obtained by performance of acts 313-315 illustrated in FIG. 3. 

FIG. 4D illustrates, in a graph, a number of excess carrier profiles for the 
doping profile of FIG. 4A, as the doping profile (see FIG. 4A) is shifted 
progressively deeper into the semiconductor. 

FIG. 4E illustrates, in a graph, a number of curves showing simulated value 
of interference measurements for diflferent intensities of the generation beam, at 
different junction depths. 

FIG. 4F illustrates, in a graph, two similar doping profiles, one identical to 
the profile shown in FIG. 4A and the other having a slightly shallower slope than 
the profile shown in FIG. 4A. 
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FIG. 5 illustrates, in a block diagram, various components used in one 
implementation of the active dopant profiler of FIG. 1 A. 

FIG. 6 illustrates, in a flow chart, another use of the measurements in 
optional act 250 (illustrated in FIG. 2 A) performed by the profiler of FIG. 1 A. 

FIG. 7A illustrates, in a graph, the active doping profile and excess carrier 
concentration as fimction of depth at various generation laser powers for a typical 
ion implant (boron, 500eV, 1x10^^ ionsW, annealed 10 seconds at 1000 »C). 

FIG.7B illustrates, in a graph, the intensity measurement of the excess 
carrier reflection as a fimction of reference arm phase for the 5, 20 and 50 mW 
power levels shown in FIG. 7A, and also the intensity measurement of the front 
surface reflectioiL 

FIG. 8A illustrates, in a prior art graph, the absorption coefficient (inverse 
of the absorption length) for silicon and amorphous silicon (a-Si) as a fimction of 
wavelength in microns. 

FIG. 8B illustrates, in a graph in accordance with the invention, the signal 
obtained in microvolts (at phase of 0°)at a generation laser power of 90 mW as a 
fimction of the thickness of the surface amorphous layer for various silicon and 
germanium ion implants. 

FIG. 9A illustrates, in a graph, the calculated excess carrier concentration 
per cm^ as a fimction of depth in microns for various values of carrier lifetime 
within an ion implanted layer approximately 0.08 microns thick. 

FIG. 9B illustrates, in a graph, the measured lock-in signal values as a 
fimction of implant dose for a variety of ion implants Usted in the key in FIG. 9C. 
In each case the set of three points indicates the signal for implant at nominal doses 
of 2.5 x 10". 5 X 10", and 1 x lO'^ ions/cm^ and implants representing doses ±5% 
above and below the nominal values. 

FIG. 9C illustrates a key for the graph illustrated in FIG. 93. 

DETAILED DFSCRIPTION 
30 A wafer fabrication system 1 00 (FIG. 1 A) in accordance with the invention 

is used to create integrated circuit (abbreviated as "IC") dice by processing a wafer 

-10- 
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(also called "semiconductor substrate") to form a "patterned wafer", measuring a 
material property of the pattemed wafer, and adjusting the processing in real time if 
necessary. The just-described processing can include annealing, and the 
measurement of a material property can be performed on a pattemed wafer after 
5 annealing, thereby to determine process conditions not obtainable by prior art 

methods, e.g. to determine anneal temperature from measurements on the annealed 
wafer. 

Measurement on pattemed wafers during fabrication as described herein 
eliminates test wafers that may be otherwise required in the prior art solely to 
1 0 monitor the fabrication process, thus reducing costs. Moreover, measurement on 
annealed wafers as described herein provides a measure of one or more properties 
that are related to the electrical characteristics (such as processing speed) of the 
devices being fabricated, because annealing results in activation of the dopants 
used in the devices. 

15 System 100 includes an active dopant profiler (also called simply 

"profiler") 103 that measures various material properties in a non destractive 
manner. Profiler 103 includes a computer 103C that is programmed with 
simulation software to generate predetermined data (see operation 1 10 in FIG. IB). 
In one implementation, operation 110 includes an act 1 1 1 wherein computer 103C 

2 0 generates a simulated value that is based on: (i) conditions present during 

stimulation, and (ii) a predetermined profile of the concentration of active dopants 
in the region under stimulation. Thereafter, in act 112, computer 103C checks if all 
predetermined profiles in a set of profiles that the wafer is likely to have been used 
in act 1 1 1 . If not, computer 103C goes to act 1 13 to change the predetermined 

2 5 profile being used, and returns to act 1 1 1 . If all profiles have been used, operation 

110 is completed. 

Wafer processing unit 101 and rapid thermal annealer 102 perform an 
operation 125 to prepare a semiconductor wafer, e.g. by forming one or more layers 
of the wafer (e.g. one of wafers 104-106). Thereafter, profiler 103 stimulates (see 

3 0 operation 126 in FIG. IB) a region in the wafer that originally has a first number of 

charge carriers, so that there are a second number of charge carriers during the 
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Stimulation (see operation 126 in FIG. IB). The stimulation can be accomplished 
in any number of ways, including e.g. by use of a beam of electromagnetic 
radiation or by a beam of electrons. 

Next, in operation 140, profiler 103 measures a property of the region that 
5 is affected by the stimulation. In one implementation, profiler 1 03 uses (see act 
141) a measurement device (such as an interferometer in one embodiment) to 
obtain one or more measured values (e.g. amplitude and phase) of a signal (such as 
a probe beam that is reflected by the charge carriers) that is affected by the 
stimulation. Next, profiler 103 compares (see act 142) the measured values with 

1 0 one or more of the simulated values (generated in operation 11 0) to identify the 

simulated value that is closest to the measured value. Thereafter, profiler 103 uses 
the predetermined profile that was used to generate the closest simulated value as 
the profile in the wafer. For example, profiler 103 determines a value of a 
semiconductor property (such as "junction depth") in the region, based on the 

15 profile used to generate the closest simulated value. Profiler 1 03 can determine the 
value in any number of ways, e.g. by computation or by looking up a table or a 
graph that relates each profile to a corresponding value of the semiconductor 
property. 

Next, in act 160, profiler 103 checks if the property value matches the 
2 0 specifications for acceptance of the wafer. If so, profiler 103 simply returns to 
operation 125, and otherwise goes to act 161. hi act 161 profiler adjusts one or 
more process conditions present during preparation of the wafer, e.g by driving a 
control signal on one or more of lines 108 and 107 to annealer 102 and wafer 
processing unit 101. Thereafter, profiler 103 returns to operation 125. 
2 5 Although m one unplementation (described above in reference to acts 111- 

1 13, the simulations are repeated for a set of profiles prior to the stimulation, in 
another implementation illustrated in FIG. IC, one or more simulations are 
repeated after the stimulation. For example, wafer processing unit 101 and rapid 
thermal annealer 102 perform operation 120, and profiler 103 performs operation 
126 (as described above in reference to FIG. IB) and thereafter performs operation 
170. In operation 170, profiler 103 uses a measurement device to obtain a 
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measured value (as described above in reference to act 141) and thereafter goes to 
act 1 72. In act 1 72, profiler 1 03 operates a simulator to generate a simulated value 
of the signal based on a predetermined profile (as described above in reference to 
act 1 1 1). Next, in act 173, profiler 103 checks if the simulated value matches the 
5 measured value (e.g. within a predetermined percentage, such as 1%). If there is no 
match, profiler 103 goes to act 174 and changes the predetermined profile that was 
used in operating the simulator (in act 172), and returns to act 172. If there is a 
match, profiler 103 goes to act 175 to determine the semiconductor property's 
value (as described above in reference to act 143). Next, profiler 103 performs acts 

10 160 and 161 as described above in reference to FIG. IB. 

In one embodiment, profiler 103 implements the stimulation operation by 
creating a number of charge carriers in a to-be-tested wafer. The charge carriers 
created by profiler 103 are in excess of a number of charge carriers (also called 
"backgroimd carriers") that are normally present in a semiconductor material (e.g. 

15 due to dopants that are defined to be atoms that occupy sites in the crystal lattice of 
the semiconductor material, and thus contribute to the electrical conductivity of the 
material) in the absence of illiunination. 

The excess carriers, when produced in one embodiment, distribute 
themselves in semiconductor material 156 (FIG. IE) in a profile 158 (see FIG. ID; 

20 defined to be the concentration in number of carriers per cubic cm) that exceeds the 
level of carriers present without stimulation (such as illimiination) formed within 
material 1 56 by the dopant atoms. Specifically, the excess carrier concentration n^ 
changes fi-om being zero outside a firont surface 153 (FIG. IE) of the 
semiconductor material 156 to a finite value inside the semiconductor material 

25 156 (thereby resulting in a step increase in the concentration at fi"ont surface 1 53), 
As the depth z ft-om fi-ont surface 153 (FIG. IE) increases, the excess carrier 
concentration increases over a depth (e.g. less than 200 A) that is at least an 
order of magnitude smaller than the wavelength (e.g. 4000 A) of probe beam 152 
within material 153 (that is the free space wavelength/refractive index). Beyond 

3 0 the just-described depth, excess carrier concentration nc changes fiirther in a 

manner proportional to a change in the concentration of dopant atoms until depth z 
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reaches the edge of the doped region, at a junction depth Zj. For example, in some 
cases, the dopant concentration rises, but in other cases the dopant concentration 
dips first and then rises, depending on the detailed shape of the doping profile. 

Beyond junction depth Zj, the excess carrier concentration returns to 
being substantially constant to a depth on the order of 10 times the depth of the 
doped region (i.e. varies less than 10%, especially for shallow doped regions less 
than or on the order of 0. 1 microns deep). The above-described profile, when 
produced as described herein, varies periodically with time, in synchronization 
with the modulation fi^uency, but does not vary periodically in space, as a 
fimction of radial distance r (FIG. ID). Instead, concentration n* simply decays 
radially (e.g. monotonically as a fimction of radial distance r) outside region 120, as 
illustrated in HG ID. Specifically, over a time period that is the inverse of the 
modulation fi^quency, profiler 103 changes concenttation ne between the values 
Uea-ncn, wherein ncn<nej<nei< Uca (FIG. ID). 

Therefore, at any given time ti, the value ne, of the carrier concentration in 
semiconductor material 156 decays as a fimction of radial distance r, without the 
creation of a wave in space (outside region 120). A profile of excess carrier 
concentration n* that is devoid of a wave (along radial distance r) is created as 
described herein when at least a majority (i.e. greater than 50%) of the charge 
carriers (defined to be both excess carriers and background carriers) that move out 
of region 120 do so due to diffusion. Within illuminated region 120 the carrier 
concentration may show spatial variation due to (1) variations in the profile of the 
generation laser beam (typically a smaller effect) and (2) doping profile variations. 

Such a "diffusive modulation" of excess charge carriers is a significant 
aspect of the invention because the predetermined data used in looking up the 
semiconductor properties are based on a diffusive solution of an equation 15 
(described below) for the movement of charge carriers fi-om region 120 (also called 
"carrier creation region"). Superposition of a wave solution (on the diffusive 
solution) degrades the accuracy of the predetermined data, because the wave 
solution perturbs the excess carrier profile away from the aperiodic profile assumed 
for the diffusive solution. 
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In one embodiment, generation beam 15 Vs intensity is modulated at a fixed 
frequency that is sufficiently low to ensure that the phase of the variation of 
concentration n^ is the same as (e.g. to within 10%) the phase of generation beam 
1 5 1 over a diffusion length (wherein diffusion length is the length over which a 
5 charge carrier decays to a value of 1/e). Therefore, concentration nc changes 

approximately synchronously with the change in intensity of generation beam 151. 
This condition ensures that the excess carrier distribution is primarily due to 
diffusion that can be modeled by a non-wave solution (rather than by a wave 
solution). 

10 To ensure the absence of a wave in space, the frequency of modulation of 

concentration n^ is selected to be several times (e.g. one or more orders of 
magnitude) smaller than the modulation frequencies used in the prior art to 
generate waves as described in, for example, U.S. Patent 4,854,710, Specifically, 
in one implementation of this invention, the modulation frequency is 

15 approximately 1 KHz that is one thousand times (three orders of magnitude) 
smaller than a 1 MHz frequency described in column 15, line 18 of U.S, Patent 
4,854,710 by Opsal. Use of such a low modulation frequency is a critical aspect in 
one embodiment of profiler 103 (FIG. 1 A), and leads to unexpected results due to 
the elimination of a wave in space, such as the "wave" described by Opsal.. 

2 0 An alternative embodiment of profiler 103 modulates a generation beam at 

a frequency that is higher than the above-discussed frequency. Specifically, the 
alternative embodiment causes periodicity in space outside of region 120. 
However, in the alternative embodiment, the wavelength of the spatial periodicity 
is selected to be greater than ten times the diffusion length, so that within region 

25 120 the effects of spatial periodicity are negligible. 

An increase in concentration n^ (as illustrated by profile 164 in FIG. IG) 
that occurs when the distribution of excess carriers forms in region 120 in response 
to generation beam 151, results in a proportional increase in the index of refraction 
"n" of the semiconductor material 1 56. An index of refraction gradient with 

30 respect to depth is thereby formed in region 120. Because of the linear 
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proportionality between excess carrier concentration and the index of refraction, 
profile 164 also represents the index of refraction as a function of depth. 

Specifically, excess carrier profile 164 (FIG. IG) may be modeled as a set 
of thin layers of excess carriers, such as layers 164A-164T A<J<T, T being the 
5 total number of layers, wherein each layer 1 64 J has a constant carrier 

concentration. The index of refraction of silicon is a function of the carrier 
concentration, so each layer 164J has a slightly different index of refraction. 
Consequently, a small amount of energy of probe beam 152 reflects from each 
mterface 165 J between layers 164J and 164J+1, with the amount of reflection 
1 0 being proportional to the carrier concentration in layer 1 64J that is located above 
interface 165 J. 

The sum of all such reflections within region 130 (that has the graded index 
of refraction) forms a component 163 that is smaller than another component 162 
of the reflection of beam 152. Specifically, component 162 is a reflection of beam 

15 152 from front surface 153, and is the sum of three subcomponents: (1) the first is 
due to the discontinuity between air and silicon. (2) the second is due to the 
sudden rise in doping profile at the surface. (3) the third is due to the sudden rise in 
the excess carrier concentration at the surface. The strongest subcomponent of 
reflection 162 from front surface 153 is the first, by several orders of magnitude as 

2 0 compared to the second and the third. 

Note that component 163 has a phase that is delayed with respect to 
component 162, because light reflected from the index gradient region 130 
propagates an additional distance Zj into the semiconductor and back out. 
Therefore, component 162 interferes with component 163, and the interference 

2 5 may be constructive, destructive, or a combination of both, dependmg upon the 

range of depths Zj over which the index gradient occurs. Thus, a signal obtained 
by summing the reflection component 162 (from front surface 153), and the 
reflection component 163 (from graded index region 130) contains within its 
amplitude and phase, information about the depth and shape of profile 164 of the 

3 0 excess carrier concentration. 
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System 100 includes a wafer processing unit 101 that performs one or more 



initial activities in act 201 (e.g. receiving wafers), and thereafter goes to operation 
210 (FIG. 2A) that is similar or identical to operation 110 (described above). Next 
(or immediately after act 201), the apparatus prepares the semiconductor wafer (in 
5 operation 220). In operation 220, unit 101 performs an act 221 to form doped 

regions, e.g. by operating an ion implanter 1011 to create, in a wafer 104 (FIG. 1 A), 
one or more regions (e.g. doped region 130 in FIG. IC) that have dopant atoms 
(e.g. boron atoms in silicon). Instead of ion implantation, any other process for 
creating doped regions, e.g. chemical vapor deposition, epitaxial deposition, 
10 evaporation, diffusion, or plasma deposition can be used in unit 101 (FIG. 1 A) to 
perform act 202. 

Thereafter, a patterned wafer 105 having one or more patterns of doped 
regions is transferred to a rapid thermal annealer 102 (FIG. 1 A) that may be 
included in system 100. Rapid thermal annealer (also called "annealer") 102 

15 performs an annealing act 222 (FIG. 2A), e.g. by heating wafer 105 (FIG. 1 A) to a 
predetermined temperature (also called "annealing temperature"), e.g. to remove 
damage that is normally caused by ion implanter 101 to the lattice structure of the 
semiconductor material in the doped regions of wafer 105. Instead of a rapid 
thermal aimealer, a furnace may be included in system 100 and used to anneal 

2 0 wafer 1 05 in act 222 (FIG. 2 A). 

Anneals (as illustrated in act 222) are typically done by heating the wafer 
rapidly with lamps (not shown) in annealer 102 (FIG. 1 A). The illumination by the 
lamps in annealer 102 may not be uniform, and the amount of heat that enters a 
patterned wafer 105 at any point may be a function of the thickness of dielectric 

2 5 layers (such as silicon dioxide or silicon nitride to be formed on surface 1 53), and 
the integrated circuit pattern therein. Specifically, the different layers (not shown) 
of doped regions in wafer 105 reflect different amounts of power, thereby causing 
variations in the amount of heating of wafer 105. Thus annealing of implanted 
wafer 105 may not be uniform, and the characteristics of a junction (formed at an 

30 interface between doped region 130 and semiconductor material 156 in FIG. IE at 
a depth Zj from surface 153) in annealed wafer 106 may vary from point-to-point. 
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Annealing in act 222 causes the dopant atoms (also called "dopants") to 
move into the lattice of the semiconductor material in a doped region 130, where 
the dopants act as donors (forming n-type material) or acceptors (forming p-type 
material). The extent to which the dopants incorporate into the lattice structure 
during act 222 is a function of the temperature at which and the time for which act 
222 is performed. The incorporation is more complete at a higher temperature or 
after a longer time. 

However, the dopants also diffuse (i.e. move) during act 222, thereby 
increasing the junction depth. The diffusion proceeds more rapidly at a higher 
temperature, and it is necessary to carefully control the annealing temperature. 
Therefore, a junction depth or profile of the concentration of dopants as a function 
of depth is measured after act 222, and the measurement is compared with 
predetermined information (e.g. a specification or infomiation obtained from 
profiles/junction depths of wafers known to be good) to determine a change (if any) 
15 to be made to the annealing process. Dynamic feedback of such to-be-made 

changes to the annealing process in real time as described herein improves the yield 
of good wafers obtained from annealmg in a manner not otherwise possible in the 
prior art. 

After annealing, wafer 106 (FIG. 1 A) is transferred from rapid thermal 
annealer 102 to profiler 103, and positioned therein (see act 224 in FIG. 2A). In an 
alternative embodunent, an active dopant profiler is integrated into a rapid thermal 
annealer and does not require positioning after completion of anneal. In one 
embodunent, profiler 103 is moved relative to wafer 106 instead of moving wafer 
104. 

Also, a non-amiealed wafer 105 can be used (moved via path 109 in FIG. 
1 A) as illustrated by branch 223 in FIG. 2A e.g. if dopant regions do not require 
annealing due to use of a method other than ion implantation, such as diffusion 
(wherein dopants are diffused into wafer 105 thermally, and are active, and there is 
no need to anneal out implant damage). Profiler 103 evaluates the efficacy of the 
3 0 dopants in a nonannealed wafer 1 05 in a manner similar to that described above for 
annealed wafer 106. 
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A Starting wafer 104 can also be used as illustrated by path 1 12 in FIG. 1 A 
and by branch 205 in FIG. 2A. Therefore, in the following description, the notation 
"104/105/106" is used to indicate that the description is equally applicable to each 
of wafers 104, 105 and 106. Similarly the notation "105/106" indicates description 
5 applicable to each of wafers 1 05 and 1 06. 

Next, after a wafer 104/105/106 is appropriately positioned (e.g. centered or 
aligned to a predetermined pattern located within the wafer), profiler 103 
stimulates a region 120 of the wafer, e.g. by creating (see operation 230 in FIG. 
2A) in a region 120 of the wafer, a number of charge carriers that are modulated at 

10 a predetermined frequency. The predetermined frequency is selected to ensure that 
a wave of the charge carriers is not created inside carrier creation region 120 during 
the act of measurement (see operation 240 in FIG. 2 A). For example, the 
predetermined frequency may be selected to be any frequency in conformance with 
the formula f^ (l/lrtx) where f is the frequency, and x is the lifetime of an 

15 excess charge carrier in the substrate. As profiler 103 does not use a "plasma 
wave" as described in U.S. Patent 4,854,710, profiler 103 is as effective in 
measuring a property of an annealed wafer 106 as in measuring a property of a non- 
annealed wafer 104/105. 

Profiler 103 (FIG, 1 A) measures a property (in operation 240 in FIG. 2 A) 

20 that is affected by charge carriers present in a doped region 130 (FIG. IC) in a 
wafer 105/106. In one implementation, the measured property is complex 
reflectance (that is, reflected portion's amplitude and phase), and profiler 103 uses 
the measurement to determine various properties (also called "semiconductor 
properties") such as junction depth, and the number of active dopants as a function 

25 of depth "Z" fi-om surface 153 of wafer 105/106. A function (called "active dopant 
profile") based on the measurement can be plotted in a graph as illustrated in FIG, 
4A described below. In other embodiments of operation 240, instead of complex 
reflectance, profiler 103 can measure other properties affected by the created 
charge carriers, such as the refractive index. 

3 0 One or more of these measurements can be used (see act 246 in FIG. 2 A) to 

lookup a material property from predetermined information (also called 
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"predetermined data") as described below in reference to FIG. 2C. Act 246 is an 
optional act, and is performed in one embodiment only after perfonnance of 
another optional operation 210. for generation of predetermined data (either 
empirically or by simulation or some combination thereof) in the form of 
5 measurements for wafers known to be good (i.e. expected measurements for wafers 
that fall within predetermined specifications for acceptance of wafers). Operation 
210 is described below in detail in reference to FIG. 2D-2F. 

One or more of these measurements may also be used (see act 260 in FIG. 
2A) by comparison against one or more predetermined limit(s) to determine if 

10 annealed wafer 106 conforms to the specification for such wafers. If wafer 106 
conforms to the specifications, wafer 106 is identified (in act 262) as being 
acceptable (e.g. by movement in the direction for further processing) and the 
conditions in wafer processing unit 101 (FIG. 1 A) and in rapid thermal annealer 
102 are left undisturbed. Thereafter, the above-described acts are repeated (as 

15 illustrated by branch 263) on another wafer or after fiirther processing on the same 
wafer. 

If a wafer 106 does not conform to the specifications, wafer 106 is 
identified (in act 26 1 ) as unacceptable (e.g. discarded) and optionally profiler 1 03 
adjusts (either automatically or under manual control) (1) the conditions (e.g. 
2 0 dosage of dopants) in unit 101 by driving a signal on a line 107 (HG. 1 A), or (2) 
the conditions (e.g. annealing temperature) in annealer 1 02 by driving a signal on 
line 108, or both. Then the above-described acts are again repeated (as illustrated 
by branch 263) on another wafer 106. 

As described below, the measurement performed by profiler 103 is non- 

2 5 destructive, is performed in a few square microns, and can be performed in a 

relatively short time (e.g. five seconds in one region or 50 seconds at 10 regions 
over a wafer). Measuring a property of annealed wafer 106 during (or immediately 
after) fabrication as described herein increases yield, as compared to an off-line 
measurement of a test wafer's properties. 

3 0 Prior to measuring a material property by performing operation 240, 

profiler 103 creates (see act 230 in FIG. 2A), in a region 120 (also called "carrier 
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creation region") of wafer 106, a concentration nc of excess carriers, and modulates 
concentration n^ (i.e. increases and decreases) as a function of time t. The excess 
carriers can be created by any method, although in one embodiment, the excess 
carriers are created by a generation beam 1 5 1 that may be a beam of 
5 electromagnetic radiation. In another embodiment, the source of excess carriers is a 
beam of electrons. 

In one embodiment, probe beam 152 is smaller in diameter than generation 
beam 151 (as illustrated in FIG. IE) due to the chromatic aberration of the focusing 
lens (e.g. lens 5 15 in FIG. 5). Moreover, probe beam 1 52 can have a longer 
10 wavelength than generation beam 15 1 , to ensure that the rate (also called 

"generation rate") of generation of carriers due to probe beam 152 is significantly 
less than the generation rate due to generation beam 1 5 1 . In one embodiment, 
generation beam 151 has a first wavelength X,g and probe beam 152 has a second 
wavelength Xp^ the second wavelength Xp being determined fi"om the formula: 



wherein ap and ag are the absorption coefficients in semiconductor material 156 
(FIG. IC) of probe beam 152 and generation beam 151 respectively, Pp and Pg are 
the powers of probe beam 152 and generation beam 151 respectively, and Wg and 
Wp are the radii of the focal spots at front surface 153 of beams 152 and 151 

2 0 respectively. This formula ensures that at front surface 1 53 the generation rate due 

to generation beam 151 is at least an order of magnitude greater than the generation 
rate due to probe beam 152. 

The wavelength of probe beam 152 is typically longer than the wavelength 
of generation beam 151 as illustrated in FIG. IE. Since, for a lens the focal spot 
25 size at surface 153 is proportional to the wavelength, probe beam 152 will typically 
focus to a larger spot size than generation beam 151. It is desirable to have the 
opposite relationship, in which probe beam 152 is smaller than generation beam 
1 5 1 at the focus, so that probe beam 1 52 is positioned within generation beam 151, 
as shown in FIG. IE. This makes the measurement less sensitive to the radial 

3 0 decay of the excess carrier concentration. 



15 



>.g > [ (10 apPp>.p)/( agPg) ] [Wg/wp]^ 
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An appropriate relationship between beams 151 and 152 is achieved using 
the chromatic aberration of a focusing lens 5 1 5 (FIG. 5). At the focus of probe 
beam 152, generation beam 151 is slightly out of focus, and, hence, has a slightly 
larger diameter than at focus. Measurements are made with this focal arrangement, 
5 m which generation beam 1 5 1 spot is at focus, with a minimum diameter, and 

generation beam 151 spot is out of focus, having a slightly larger diameter than the 
probe beam 152. 

In one embodiment, profiler 103 implements the above-described act 230 
(FIG. 2A) by: generating (act 23 1 in FIG. 2B) a beam 151 (FIG. IE) of photons 
10 that have energy greater than the bandgap energy of the semiconductor material in 
doped region 130, modulatmg (act 232 in FIG. 2B) beam 151 at a frequency 
selected to avoid the creation of a wave (as described above), and focusing (act 233 
in HG. 2B) beam 151 on doped region 130. However, in an alternative 
embodiment described below in reference to measurement of the depth of an 
15 amorphous layer, profiler 103 implements the above-described act 231 (FIG. 2B) 
by generating a beam 1 5 1 of photons that have energy lower than the bandgap 
energy of the semiconductor material (amorphous silicon) in doped region 130. 

Dependmg on the implementation, profiler 103 modulates the intensity of 
generation beam 1 5 1 at any frequency in the range of 1 Hz to 20,000 Hz. The 
20 modulation fi-equency can be, for example, 1000 Hz, and may require at least 10 
cycles for a lock-in amplifier to generate a reflectance measurement (based on a 
probe beam as described below in reference to act 242), or 10 milliseconds to 
perform each reflectance measurement. In one example, the throughput is 30 
wafers per hour, or 120 seconds per wafer, with each wafer having a measurement 
25 taken in at least ten regions. 

If a material property measurement requires several reflectance 
measurements (e.g. a single region 120 requires a number of reflectance 
measurements for each of a corresponding number of average carrier 
concentrations that may be obtained at a range of powers of generation beam 1 5 1 
(10 powers are linear spaced between 5 mW and 100 mW average power, for 
example)), profiler 103 takes several seconds (e.g. 10-100 seconds) for each wafer 
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104/105/106. Hence, the 10 millisecond speed of reflectance measurement per 
region allows for real time control in the fabrication of wafers by apparatus 1 00 
(FIG. 1 A) using method 200 (FIG. 2 A). 

In another implementation of act 230, instead of using beam 1 5 1 of 
5 photons, profiler 103 uses a beam of charged particles, such as electrons or ions. 
The beam of charged particles is modulated and focused in the same manner as that 
described herein in reference to beam 151 to generate the charge carriers in doped 
region 130. Instead of a beam of photons or a beam of electrons, any other 
mechanism (such as a combination of photons and electrons) can be used to create 

10 charge carriers in act 230 (FIG. 2A). 

In act 240, one implementation of profiler 103 focuses (see act 242 in FIG. 
2A) on a region (also called "carrier creation region") 120 illuminated by beam 
151, another beam 152 (FIG. IE) that is used to detect the number of charge 
carriers in wafer 104/105/106 when illuminated by beam 15 1 . In one embodiment, 

15 beam 152 (also called "probe beam") contains photons having energy lower than 
the bandgap energy of the semiconductor material in carrier creation region 120, 
Such a probe beam 152 limits the creation of additional carriers (due to the probe 
beam, also called "measurement-related carriers") when beam 152 is incident on 
carrier creation region 120, thereby to maintain the charge carrier concentration 

20 approximately the same prior to and during measurement (see act 243 in FIG. 2A) 
of an interference signal as described below. 

Next, profiler 103 measures (see act 243 in FIG, 2A) the amplitude and 
phase of a signal generated by interference between a probe beam 152 (FIG. IE) 
reflected by the excess charge carriers within region 156 (FIG. IE) and a reference 

2 5 beam (that may be either a portion of the probe beam reflected fi-om front surface 

1 53, or another portion of the probe beam that has a variable phase, the phase being 
varied as described below in reference to act 244). As the interference signal being 
measured is modulated at the frequency of modulation of the charge carriers in 
carrier creation region 120, a lock-in amplifier (described below) may be used to 

3 0 improve accuracy of the measurement. 
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The measurement in act 243 (FIG. 2A) provides an indication of an average 
concentration nav of cliarge carriers in doped region 130 near surface 153, wherein 
the average concentration nav is a root mean square average that is measured over 
the period of one (or more) modulation cycle(s) at the modulation frequency of 
5 generation beam 151. Concentration nav in turn indicates, under certain conditions 
as discussed below, a material property, e.g. the junction depth in doped region 
130. 

In one embodiment, the location at which the charge carriers are created is 
not changed between two or more measurements. Instead, in one implementation, 
10 profiler 103 performs a number of measurements at the same location (e.g. at least 
two measurements for two different powers of generation beam 151) in wafer 
105/106, but changes a parameter used to create the charge carriers. The parameter 
can be, for example, the average carrier concentration nav in region 120. 

Concentration nav is changed e.g. by changing the intensity of generation 
15 beam 1 5 1 (e.g. by changing the power or the diameter), and act 243 is repeated. 
Alternatively, profiler 103 can change the location of carrier creation region 120 
and perform a number of such measurements. In one implementation, the locations 
of each of probe beam 152 and generation beam 151 are changed to obtain a linear 
scan across a wafer 104/105/106, while holding the beams 151 and 152 coincident 
each with the other. Also, instead of or in addition to act 241, profiler 1 03 changes 
a parameter used in the measurement as illustrated by act 244 in FIG. 2A, e.g. 
changes phase of the reference beam. 

In the above-described embodiments, a probe beam 152 having photons of 
energy below the bandgap energy of wafer 156 is used, although in another 
embodiment probe beam 152 has photons of energy equal to or sUghtly above (e.g. 
5% above) the bandgap energy. Certain additional carriers (called "measurement- 
related carriers") created by probe beam 1 52 are in a sufficiently small percentage 
(e.g. an order of magnitude smaller than the number created by the generating 
beam) to provide a reasonably accurate measurement (e.g. to within 5%). Note that 
the overall accuracy of a measurement as described herein is also governed by 
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Other inaccuracies involved in the act of measuring, e.g. inaccuracies in a 
measurement device, such as a lock-in amplifier. 

Therefore, in one embodiment the inaccuracy caused by the measurement- 
related carriers is kept only as small as necessary to maintain the overall accuracy 
5 below a predetermined limit. Specifically, the percentage of measurement-related 
carriers is kept sufficiently small when the rate per unit volvune of the carriers 
generated by generation beam 151 (obtained by dividing the photon flux per unit 
area by the absorption length), is at least one order of magnitude (or more) larger 
than for probe beam 152. 

10 The photon flux per unit area in region 120 due to generation beam 151 is 

the number of photons per unit area obtained by dividing the power P of generation 
beam 151 by the area (tiWq^) of illumination, where Wo is the radius of generation 
beam 1 5 1 , by Plank's constant h and the ratio of the speed c of light to the 
wavelength X as shown in the following formula: photon flux = (P/7tWo^)x 

15 (l/h(c/X)). The absorption length is the depth from surface 153 at which the 

intensity of generation beam 151 drops to (1/e) of the intensity at surface 153 (see 
equation 23). 

In one implementation, probe beam 152 has a generation rate one or more 
orders of magnitude smaller than the generation rate of generation beam 151. As 
2 0 noted above, the difference in generation rates is obtained by using beams 151 and 
152 that have different absorption lengths in the semiconductor material of wafer 
156, or by generating beams 151 and 152 at different powers or different diameters, 
or all of the above. In various implementations, the pair of beams 151 and 152 are 
generated by one of the following pairs of lasers: (AlGaAs, InGaAs), (Ar, InGaAs), 

2 5 (Nd: YAG, hiGaAs), and (Nd: YAG, AlGaAs). 

In one or more of the implementations, e.g. for use of lasers (NdiYAG, 
AlGaAs), the power of probe beam's laser (e.g. AlGaAs) is maintained less than the 
power of generation beam's laser (e.g. Nd:YAG) because the absorption length of 
the probe beam is a firaction (e.g. one-tenth) of the absorption length of the 

3 0 generation beam. In another example, a probe beam 1 52 formed by a HeNe laser is 

maintained at a power less than or equal to 1/4^** power of generation beam 151 
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formed by an Ar laser (having an absorption length 1 .2 ^m that is 1/4*^ the S.Oum 
length of the HeNe laser beam). In the just-described implementation, the power 
of the reflected portion of probe beam 152 is maintained large enough (by having a 
sufficiently large power of probe beam 152) to be detected with sufficient accuracy 
5 (e.g. with error of 5% or less) required for reflectance measurements as described 
herein. 

In one variant of this implementation, the difference between the generation 
rates of beams 151 and 152 is one order of magnitude only at surface 153 (FIG. 
IC). In a second variant, the order of magnitude difference is maintained 

10 throughout junction depth "Zj" of doped region 130 in wafer 105/106. e.g. 
throughout depth of 0.3 microns. In a third variant, the order of magnitude 
difference is maintained throughout a predetermined fraction (e.g. 1/2) of the 
junction depth Zj. 

In one embodiment, each measurement for a wafer 104/105/106 is 

15 compared (in act 260 in FIG. 2A) with a predetermined range, and if any 

measurement falls outside the range, the wafer is rejected, hi one implementation, 
computer 103C displays on monitor 103M a message indicating that measurements 
identify a wafer 104/105/106 as unacceptable, while in another implementation 
computer 103C drives a signal to a robot (not shown) to move wafer 104/105/106 

2 0 into a bin of rejected wafers (if rejected). The acceptable wafers are processed 
further in the noimal manner (see act 262 in FIG. 2A). 

In addition, act 240 is used in one implementation to screen out starting 
wafers formed of bare silicon. When defects in such bare silicon are identified at 
the beginning, the method results in correction of the wafer fabrication process to 

2 5 ensure a sufficiently low defect level and eliminate the cost and use of a starting 

wafer 106 formed of epitaxial material. Starting wafers formed of pure sihcon 
(also called "prime wafers") are processed by profiler 103 in a manner identical to 
starting wafer 104 as described herein. 

Two or more of the interference measurements made in act 243 (FIG. 2A) 

3 0 may be used (in an operation 250) to look up a material property of wafer 

104/105/106. For example, in the front surface embodiment, profiler 103 
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performs acts 251-253 illustrated in' FIG. 2C. Specifically, in act 251, profiler 103 
transforms the interference measurement, e.g. by multiplying with a constant. Note 
that act 251 is an optional act, for example if the predetermined data is scaled. 

Next, in act 252, profiler 103 performs look up of predetermined data by 
5 determining a location 2551 (FIG. 2D) of the measurement (either the raw 

measurement or the scaled measurement, depending on whether the predetermined 
data is scaled) on a line 255. In performing the lookup, the measured amplitude is 
used with a sign that is determined from the measured phase, e.g. if the phase 
measurement is less than 180 degrees, the measured amplitude is used as a positive 

1 0 value, and otherwise as a negative value. Note that a specific value of the phase 
measurement is not used in this embodiment (other than to determine the sign), 
although such a value is actually used in another embodiment (described below in 
reference to FIG. 6). 

Profiler 103 uses a line (also called "curve") 255 that is a plot (along the 

15 y axis) of the measured signal (e.g. in microvolts) as a function of a material 

property (along the x axis), such as junction depth (for a selected concentration nc 
of excess carriers, as determined in this example, by generation beam power of 
lOmW). Specifically, in act 253, profiler 103 reads ofiFa value vl (also called 
"first value") of the material property from the determined location. In the 

20 example illustrated in FIG. 2D, signal SI (obtained by interference between the 

reflections by front surface and excess carriers) has a value of 10 |iV, and profiler 
1 03 finds a value vl = 0.04 [xm for the junction depth. 

In the just-described example, programmed computer 103C can compare 
either SI or vl (in act 260 illustrated in FIG. 2 A) with an appropriate one of ranges 

25 40 to -20 and 0.03 to 0.05 |am to determine that the wafer is acceptable (in act 262). 
Profiler 103 (FIG. 1 A) may perform act 241 (by increasing the generation beam 
power to 40 mW), and repeat acts 252 and 253, to fmd a second value v2 of the 
material property from another line 256 (FIG. 2E). In the example illustrated in 
FIG. 2E, profiler 103 uses an interference signal S2 of value 150 |iV to again find 

30 a value v2 = 0.04 |im for the junction depth. Also, in the example in FIG. 2E, 
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signal S2 is obtained by increasing the intensity of generation beam 151 to 40 mW 
(from the 1 0 mW used to obtain the interference signal of value SI). 

Graphs 2D-2F are computed using one or more methods of the type 
described below in reference to FIG. 3, wherein a doping profile is assumed. 
5 Simulated values for both amplitude and phase (e.g. phase of 0° and 1 80° results in 
positive and negative values respectively for the intensity) are then calculated (e.g. 
see act 3 15 described below) as a function of generation laser power and profile 
depth to obtain graphs 2D-2F. In practice, the real doping profile may deviate from 
the assumed profile. In this case, values vl , v2 and v3 may not be identical, and 
10 the three values may be averaged. 

In one embodiment, the result is reported as the average va of the three 
values vl-v3, and the range vr (obtained as the difference between the maximum 
and the minimum among values vl-v3). In another embodiment, the result is 
reported as the average va of the three values vl-v3 and a standard deviation vsd 
15 from the average va. That is, vsd = square root (((va-vl)^ + (va-v2)^ + (va-v3)^ 

)/9). The range or the standard deviation are also compared with acceptable values 
for range or standard deviation as may be provided in a manufacturing 
specification. 

Measuring junction depths as described above provides an unexpected 
2 0 result, considering that at least one prior art reference, namely U.S. Patent 

4,854,71 0 granted to Qpsal teaches that depth information cannot be obtained in 
the absence of a plasma wave (specifically, Opsal states in column 4, lines 33-35, 
"[hjowever, in applications where sample variations as a function of depth need to 
be studied, it is necessary to generate and study plasma waves"). 

Graphs (e.g. see lines 255-257 m FlGs. 2D-2F), that are used to determine a 
material property (or a process condition) can be generated in any way. In a first 
embodiment, a set of wafers (also called "reference wafers") is selected or prepared 
to have a range of material properties (by varying process conditions, such as 
implant energy, dose or anneal temperature), and thereafter profiler 103 is used to 
obtain interference measurements as described herein, and generate best-fit lines 
for each of the measurement conditions (e.g. for each generation beam power as 
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described above). In a second embodiment, a number of wafers (also called 
"reference wafers") are subjected to intensit>' measurements in profiler 103 (as 
described above), followed by use of a conventional measurement technique, such 
as spreading resistance profiling (abbreviated as "SRP") to determine the actual 
5 doping profile therein. 

In both embodiments, programmed computer 103C generates each of lines 
255-257 (FIGs. 2D-2F) from intensity measurements under the same conditions 
(e.g. diffiisive modulation) on the set of reference wafers having different material 
properties (e.g. junction depths in the range of 0.02 [im to 0,08 fim, at increments 

10 of 0.005 |im). The reference wafers may be prepared by ion implantation of 

species such as Boron, Arsenic, Phosphorous, or BF2 at an energy range of 0.2 to 
5KeV and dosage of ixlO^Vcm^, followed by annealing (e.g. for 10 sieconds) at 
each of the temperatures in the range 900-1 050°C, with 50 ®C increments. 

Thereafter, if the material properties are not known, SRPs are prepared by 

15 breaking the wafers to expose the ion-implanted layer followed by beveled lapping 
and probing to measure the profile of the concentration of active dopants as a 
fimction of depth. Therefore, at the end of the preparation of SRP, the graphs (e.g. 
FIG. 4A) provide a plot of the active dopant concentration (atoms/cm'^) along the y 
axis as a function of depth (in microns) along the x axis. 

20 Therefore, profiler 103 obtains a number of measurements SA-SN (FIG. 

2D) for each of points 255A-255N (A<I<N, N being the total number of 
measurements), with beams 151 and 152 coincident in the same region 120 (FIG. 
IC) on each of the reference wafers. Thereafter, profiler 103 fits points 25 5 A- 
255N to a curve 255 (e.g. represented by a linear approximation of the form y = - 

25 3.016 X + 130.88 in FIG. 2D). In a similar manner, profiler 103 generates points 
256A-256N and 257A-257N for generation beam powers of 40 mW and 100 mW 
respectively, and thereafter generates the corresponding curves 256 and 257. Note 
that each of curves 255-257 is a sinusoidal curve that can be approximated by a 
straight line (as shown in the example in FIG. 2D) over a portion of the curve. 

3 0 Alternatively, sinusoidal curves 255-257 that are obtained fi-om curve 

fitting of points 255A-255N, 256A-256N and 257A-257N may be used directly. 
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without a linear fit to determine a straight line. Use of sinusoidal curves 255-257 is 
more accurate than use of a linear approximation, although the linear 
approximation is simpler to implement. Also, instead of approximating a curve 255 
with a straight line equation, a second or higher order differential equation can be 
5 fitted to the points, and the differential equation may be used to obtain the property 
measurement (in the manner described herein, as would be apparent to a person 
skilled in the art of computer programming). Moreover, instead of a reflectance 
measurement being used to measure a property of the semiconductor material, a 
change in the index of refi-action can also be used in a similar manner. 

^0 After one or more of the above-described graphs (see FIGs. 2D-2F) are 

prepared, the material properties of a wafer under fabrication are determined by the 
above-described method 200 (FIG. 2A) without the need to break and lap the 
wafer, because profiler 103 simply uses the above-described graphs to generate 
measurements of material properties. Therefore, profiler 103 eliminates the cost 

15 associated with test wafers otherwise required by the prior art methods (for 
breaking and lapping). 

Although in the above description, computer 1 03C has been described as 
performing various computations for the preparation of curves (e.g. curve 255 in 
FIG. 2D) used to measure material properties, such graphs can be prepared by 

2 0 another computer, or alternatively can be prepared by manually performing the 
above-described acts. 

Moreover, although in one embodiment the above-described curves (e.g. 
FIGs. 2D-2F) are drawn, in another embodiment such graphs are not prepared and 
instead the reflectance measurements are simply used to perform the various acts of 

2 5 method 200, by use of equations related to such graphs. For example, instead of 

drawing a curve 255 (FIG. 2D), an equation for the curve is determined by fitting 
(as discussed above), and thereafter the equation is used to obtain the material 
property. 

hi one implementation, probe beam 152 (FIG. IE) is a laser beam having a 

3 0 wavelength greater than 1 .05 (the wavelength at which photons have 

approximately the same energy as the bandgap energy of silicon). Note that the 
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wavelength of probe beam 152 depends on the bandgap energy and therefore on the 
specific material in wafer 105/106, and is different for germanium. 

In one embodiment, the predetemiined data for use in act 246 (described 
above) is prepared by performing method 3 1 0 (FIG. 3) that implements act 1 11 . 
5 Specifically, in act 3 1 1 , computer 103C (FIG. 1 A) receives (see act 3 1 1 in FIG. 3) 
a dopant profile 401 (FIG. 4A) for a selected junction depth (e.g. a depth of 0.035 
|im illustrated in FIG. 4 A). Dopant profile 401 may be obtained by either 
simulation (e.g. by use of the simulators Atlas and Athena, both available from 
Silvaco International, Santa Clara, California) or by a conventional method such as 

1 0 spreading resistance profile (SRP) or secondary ion mass spectrometry (SIMS). 

Next, computer 103C uses (see act 312 in FIG. 3) a simulator to determine a profile 
402J (FIG. 4B) of excess carriers for a selected average concentration of excess 
carriers (e.g. as determined by a selected power of a generation beam), e.g. 10 mW. 
Thereafter, for each of a number of points 403JA-403 JR (wherein A<K<R, 

15 R being the total number of such points) computer 103C multiplies (as illustrated 
by act 3 13 in FIG. 3) the foUov^dng two multiplicands to obtain a product: (1) 
derivative of profile 402J with respect to depth from front surface (i.e. dy/dz) and 
(2) cos (2knz), wherein z is the depth, kr=2n/X, X is the wavelength of the probe 
beam, and n is the index of refraction of the substrate. 

2 0 Next, computer 1 03C integrates (see act 3 1 4 in FIG. 3) the products (at 

each of points 403JA-403JR) with respect to depth z, adds the value of the excess 
carrier concentration at the surface, and multiplies (see act 315 in FIG. 3) the result 
of integrating with a constant. The constant is based on calibration of profiler 103, 
and provides the simulated value of the signal obtained by interference between 

2 5 portions of probe beam 1 52 that are reflected by the excess carriers and by front 

surface 1 53. Computer 103C may (in an optional act not illustrated in FIG. 3) plot 
a point 404J (FIG. 4C) in a graph that shows the simulated value along the y axis as 
a function of generation beam power along the x axis. 

Then, computer 103C checks (in act 316 in FIG. 3) if a number of selected 

3 0 powers of generation beam 1 5 1 (FIG. 1 E) have been processed, as described above 

for profile 403 J. If not, computer 1 03C returns to act 3 1 2 (and performs acts 3 1 3- 
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3 1 5) for another profile (for another selected power of generation beam 151). In 
one example, computer 103C performs acts 312-315 for each power in the set of 
0.1, 0.2, 1, 2, 4, 6, 8, 10, 14, 20, 40, 60, 80 and 100 mW (a range over which 
generation beam power is expected to be varied), thereby to obtain each of a 
5 corresponding number of points 404A-404P, wherein A<J<P, P being the total 
nimiber of such profiles. 

Thereafter, computer 103C draws a line 404 (FIG. 4C) that fits the points 
404A-404P. Computer 103C repeats method 310 (i.e. acts 311-316 described 
above) for each of a number of different junction depths. Note that instead of 
1 0 generating depth profiles for each of the different junction depths by SRP or by 
simulation, the excess carrier profiles of act 312 can be generated by simply 
offsetting each of profiles 402A-402P (FIG. 4B) by the increment in junction depth 
(e.g. shifting each profile 402J to the right by a distance of 0.005 ^m on the x axis, 
as illustrated by profile 403 J in FIG. 4D). Thereafter, the shifted profiles are 
15 processed as described above in reference to acts 313-315 (FIG. 3) to obtain 

another line (e.g. line 405 illustrated in FIG. 4E) that relates the generation beam 
power to the simulated value of the interference signal. In this manner several 
lines 404-409 are obtained (see FIG. 4E). 

Next, for a given power of generation beam 151 (e.g. power 90 mW), 
2 0 computer 103C reads off a number of points on lines 404-409 to determine the 

simulated value of the interference signal for each of a number of junction depths. 
For example, simulated values for points 404L-409L are read off for each of 
junction depths 0.020 - 0.080 ^im at increments of 0.005 [im. Thereafter, computer 
103C plots points 404L-409L on a graph of simulated value of the interference 
2 5 signal as a function of junction depth (e.g. see FIG. 2E, wherem points 256A-256N 
correspond to points 404L-409L), and draws a line (e.g. line 256) that best fits the 
pomts. Computer 103C prepares simUar graphs for each concentration of the 
excess carriers (e.g. for each of a number of generation beam powers), and 
thereafter uses the graphs as described above in reference to act 246 (FIG. 2A). 

Computer 103C also prepares similar graphs for doping profiles that are 
similar to original profile 401 (FIG. 4A), but are modified by making the original 
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profile deeper (e.g. see profile 410 in FIG. 4F) or shallower. In one example, an 
original profile 401 (FIG. 4 A) is again used to prepare the similar profiles that are 
deeper or shallower, in accordance with the following procedure. The depth values 
were multiplied by a constant to expand the depth along x axis (in the example, the 
5 constant was 1 .25). This makes the scaled profile (not shown) deeper, since each 
doping point is L25 times deeper than it was in the original profile 401. The depth 
values are then shifted by a constant (e.g. in this case 70 angstroms) so that the 
modified profile 410 (FIG. 4F) begins at the same point as original profile 401 . 
In one embodiment, a line 255 (FIG. 2D) showing interference signal vs. 

10 generation beam power is obtained, e.g. by measurements at the center of a 200 
mm diameter wafer. The interference signal is then measured at a number of 
points, e.g. at points that are spaced 1 mm apart, at each of three generation beam 
powers of 1 0, 40 and 1 00 mW, along a diameter scan of the wafer. Lines 255-257 
(FIGs. 2D-2F) are thereafter used to determine the doping profile within the wafer 

15 at each point within the scan. In this manner, the profile across the wafer may be 
determined to be uniform (e.g. falls within certain preset limits, such as 10% 
variation in the depth at v^ich a concentration of 10^* dopant atoms per cubic 
centimeter is present). 

Moreover, shifts in peak doping level are equivalent to a change in the 

20 slope of the profile 401 (FIG. 4 A) of active dopants. For example, a 10% increase 
in peak doping is equivalent to a 10% increase in the slope of profile 401 (FIG. 
4A). Hence, the peak doping sensitivity of the interference signal can be 
characterized, and used to measure the peak doping in a manner equivalent to that 
described above in reference to FIG. 3 for measuring a change in the junction 

25 depth. 

In the first embodiment wherein a portion of probe beam 1 52 reflected by 
front surface 153 (FIG. IE) interferes with another portion reflected by excess 
carriers, probe beam 152 is generated by a laser 501 (FIG. 5), that can be a 
conventional laser diode, such as a 1 .48 mm wavelength InGaAs diode with a 
3 0 maximum power of 70 m W made by Hewlett-Packard of Palo Alto, CA. 
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In a second embodiment wherein probe beam 152 is interfered with a phase 
variable beam, laser 501 is a distributed Bragg reflector (DBR) AlGaAs laser with 
a wavelength of 1083 nm and a power of 50 mW (Spectra Diode Labs, San Jose, 
CA). 

5 A DBR laser is used in the second embodiment because it has a coherence 

length in excess of a meter. This simplifies interferometer design, since the 
reference beam length is not critical as long as the difference in path length 
between the reference beam path and the probe beam path is shorter than the 
coherence length (the probe beam path length is twice the distance from beam 

10 splitter 512 to the wafer 516; the reference beam path length istwice the distance 
from beam splitter 512 to mirror 513). 

The output of laser 501 is collimated using lens 502 to provide a coUimated 
beam 503 with a diameter of 3 mm. Lens 502 can be, for example, part number 
WT-CY3-163-10B-0.5 available from Wave Optics, Mountain View, California. 

15 In one embodiment, a generation beam 1 5 1 is created by an above bandgap laser 

505, such as an AlGaAs diode laser v^th a wavelength of 830 nm and power of 200 
mW, available from Spectra Diode Labs, San Jose, California. Profiler 103 
includes a lens 507, which is part number 06GLC002/810 available from Milles 
Griot Corporation, Irvine, California. Lens 507 coUimates the beam from laser 505 

20 to generate a collimated beam 151 with a diameter of 3 mm. Lens 507 is mounted 
on a positioner (not shown) for providing motion to beam 151 with respect to beam 
152. 

The relation between wavelengths of beams 151 and 152 produced by lasers 
501 and 505 is a critical aspect in one embodiment and leads to unexpected results, 
25 for example when beam 151 contains photons having energy above silicon's 

bandgap energy and beam 152 contains photons having energy approximately the 
same as or less than the bandgap energy. In this example, for a silicon wafer the 
830 nm and 1083 nm wavelength beams provide one or more benefits described 
herein. 

3 0 Wavelength 830 nm is considered particularly suitable for generation beam 

151 because the absorption length in silicon is about 15 microns. Thus, the 
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absorption length is much greater than the junction depth, and creation of excess 
charge carriers is nearly uniform over the depth of concern in the measurement. 
Because the photon energy is close to the bandgap energy, photon generation is 
more efficient, with less energy going directly into heating the semiconductor. 
5 Also, the absorption length at wavelength 1083 nm is about 300 microns, 

and therefore the number of excess carriers being created by such a probe beam is 
sufficiently low to ensure minimum perturbation to the excess carrier distribution. 
Moreover, the absorption length at wavelength 1083 is short enough that very little 
reflection firom a back surface of the wafer is seen (wafers are typically 600-800 
10 microns thick), since the back surface reflection can potentially cause spurious 
signals. 

Beams 151 and 152 are combined using dichroic splitter 510 (such as a 
partially transmissive mirror (e.g. part number 1918-b available from Dominar of 
Santa Clarzi, CA), forming a superposed beam 511. Beam 511 passes through a 

15 50:50 beam splitter 512 (e.g. part number 2005 from Dominar ) that directs a 

portion of beam 51 1 to detectors 522a and 522b (via filter 520 and polarizing beam 
splitter 521), for use in measurement of an interference signal. The remainder of 
beam 51 1 passes through a 90:10 beam splitter 514 (available from Precision 
Applied Products of FuUerton, CA, by specifying 93.3% transmission at 0.83 

2 0 microns wavelength and 90% transmission at 1 ,48 microns wavelength), and an 
objective lens 515 (such as a lOOX, 0.8 NA lens made by Olympus of Tokyo 
Japan). Objective lens 515 focuses the combined beam 511 onto wafer 516. 

Note that the specifications for beam splitter 514 are selected based on the 
wavelengths of the generation and probe beams to ensure that a majority of the 

2 5 power is transmitted and a smaller amount (e.g. 1 0%) of the power is reflected. 

Note also that probe beam 152 is focused only in a carrier creation region 120 
(FIG. IE) that is formed by focusing generation beam 151 (FIG. 5). Specifically, 
because of chromatic aberration, the focal planes of beams 151 and 152 differ 
slightly. The size of the focal spot for beam 152 is smaller than the size of the 

3 0 focal spot for beam 1 5 1 by virtue of the shorter wavelength of beam 151. If wafer 

516 is placed in the focal plane of beam 152, beam 151 will be slightly out of focus 
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and its spot on front surface 153 (FIG. lE),of wafer 516 (FIG. 5) will be larger in 

diameter and fully overlay the focal spot of beam 152. 

Light reflected from wafer 516 passes back through objective lens 515, 

90:10 beam splitter 514, and into 50:50 beam splitter 512. Half of the light 
5 reaching beam splitter 5 12 is directed back through filter 520 (which is a bandpass 

filter that blocks the light from beam 1 5 1 but passes the light from beam 1 52). 

Filter 520 can be, for example, Schott glass RG830, available from Spindler & 

Hoyer Corporation of Goettingen, Germany. Alternately, filter 520 can be a 

narrow-band pass filter with a center wavelength of 1080 nm, available fiom 
10 MellesGriotof Irvine, CA. 

Filter 520 removes photons of generation beam 151 from the reflected 

beam, thereby allowing detector 522a to see only the photons of probe beam 1 52. 

Filter 520 is a critical component in one embodiment and provides the unexpected 

result of eliminating feed-through of the modulated signal (generated by beam 151) 
15 to detector 522a that would otherwise be present when using a prior art system. In 

this particular implementation, germanium is used in photo detector 522a to 

provide sensitivity to photons of wavelength 1083 nanometers that are generated by 

laser 501. 

In the first embodiment, a reference beam is formed by a portion of probe 
2 0 beam 1 52 that reflects fix>m front surface 1 53 (FIG. 1 E), and 50:50 beam splitter 

5 12 diverts 50% of the reflected beam from fi-ont surface 1 53 toward detector 522a. 

Note that in the first embodiment, beam splitter 521, detector 522b, and amplifier 

523b are not used (i.e. are not present). 

Detector 522a is a photocell (such as a photodiode or a phototransistor, e.g. 
25 J 1 6-8SP-R05M-HS fi-om EG&G Judson of MontgomeryviUe, PA, USA) that 

converts the incident interference signal into a current. Amplifier 523a converts 

the current to an amplified current which is then sent to an amplifier 524 that in 

turn is coupled to a lock-in amplifier 525 (such as model 830 available from 

Stanford Research Systems, Suimyvale, California). 
! 0 Lock-in amplifier 525 includes a reference oscillator at the lock-in detection 

frequency. This oscillator is coupled to a laser driver 526 to provide a signal to 
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laser 505 that is modulated at the same frequency as the signal provided by lock-in 
amplifier 525. Lock-in amplifier 525 provides a signal indicating the amplitude as 
well as phase of reflected beam with respect to modulation by laser driver 526 to a 
processor 527, such as a personal computer running software to capture and display 
5 the signal in an appropriate manner (e.g. in a graph). The signal may also be stored 
in the personal computer (e.g. in a database on the hard disk) for later processing. 

In one implementation, personal computer 527 has a line 528 that is 
coupled to lines 107 and 108 (described above in reference to FIG 1 A) thereby to 
control the acts performed by ion implanter 101 and rapid thermal annealer 102 
10 based on measurement of one or more material properties as described herein. 

Beam splitter 514 diverts 10% of the return beam from wafer 5 16 via a lens 

517 (such as tube lens 81845 available from - Nikon of Tokyo, Japan) to a camera 

5 1 8 (such as a CCD camera, e.g. model 85400 available from FJW Industries of . 
Palatine, XL). The signal provided by camera 5 1 8 is fed into a vision system (not 

15 shown in FIG. 5), such as model ASP-60CR-1 1-S available from Cognex 
Corporation, Boston, Massachusetts. 

Positioning of wafer 516 with respect to the combined beam 51 1 is 
accomplished using a microscope that includes stage 529, objective lens 515, beam 
splitter 514, lens 517 and camera 518. Stage 529 can be used is used to move 

2 0 wafer 1 06 relative to beam 511 in the X, Y and Z directions. Specifically, stage 

529 can be used to move wafer 516 in the vertical direction along the Z axis to 
adjust focus, and in a horizontal plane to adjust the position of region 120 of FIG 
IE relative to beam 511. 

In a second embodiment, an independent beam 531 having a variable phase 
25 is used as a reference beam (instead of using just portion 1 62 that is reflected by 
front surface 153). Specifically, reference beam 53 1 is a portion of probe beam 
152, and compensator 506 rotates the polarization 90°, so that the polarization of 
reference beam 531 is orthogonal to the polarization of portion 162 reflected from 
wafer surface 153. In one implementation, compensator 506 is model 5540, 

3 0 available from New Focus Inc., Santa Clara, California. The phase of reference 

beam 531 is varied by adjusting the path length using piezoelectric positioner 504 
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that is located behind mirror 513. The phase and polarization of beam 531 can be 
changed independent of probe beam portion 162. 

A measurement (also called "front surface reference") of interference 
between the reference beam and reflection from the front surface is a non-wave 
signal (also referred to as a "dc» signal) that varies at the same rate as variation of 
the length of the path of the reference beam (also called "reference arm length"), 
e.g. by a piezoelectric device 504 that moves mirror 513. In the embodiment 
illustrated in FIG. 5, the front surface reference measurement is provided by sum 
detector 522a. 

Another measurement (also called "excess carrier reference measurement") 
of interference between the reference beam and reflection by the excess carriers is 
signal that is modulated at the same frequency as generation beam 1 5 1 and that is 
delayed in phase as compared to the front surface reference. The- difference 
between the two measurements provides the phase difference that indicates the 
absolute junction depth as described below, hi the embodiment illustrated in FIG. 
5, the excess carrier reference measurement is provided by lock-in amplifier 525. 

Specifically, in the second embodiment, the path length (also called 
"reference arm length") of reference beam 531 (to and from front surface 153) is 
set to be an integral multiple of the probe beam wavelength. This is done by 
adjusting the position of mirror 513 by applymg a voltage to piezoelectric element 
504 so that the signal in sum detector 522a is at maximum. The maximum occurs 
when there is constructive interference between beam 531 and reflection 162 of 
probe beam 152 from front surface 153. 

Next, measurements at sum detector 522a and at lock-in amplifier 525 are 
2 5 recorded (note that the lock-in amplifier signal corresponds to the difference 

between the outputs of detectors 522a and 522b). Thereafter, the reference arm 
length is incremented a small firaction of a wavelength, and these two 
measurements are again recorded. The reference arm length is changed by 
changing the voltage applied to piezoelectric element 504 that controls the location 
30 ofmirror513. 



20 
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Such measurements may be repeated over several intervals of wavelength 
of probe beam 152, so that a smooth curve is obtained (over a half-wavelength 
interval; a curve having greater than 10: 1 signal-to-noise ratio is considered smooth 
in one embodiment). The measurements by each of lock-in amplifier 525 and 
5 detector 522a are plotted to obtain two curves, one for each of the front surface 
reference measurement and the excess carrier reference measurement that are 
plotted against normalized reference arm length (obtained by dividing the reference 
arm length by the wavelength of probe beam and multiplying by the refractive 
index of the semiconductor material) as illustrated in FIG. 7B. 

10 Thereafter, the phase difference G is measured, e.g. to be 73 degrees, and 

the phase difference when divided by 360 degrees and multiplied by a constant 
(probe beam wavelength divided by index of refraction in the semiconductor) 
yields the absolute value of the junction depth (the depth at which the 
concentration of dopants is equal to a predetermined concentration, e.g. 10* Vcm^). 

15 The phase shift of 73° corresponds to a jimction depth of 0.044 ^m, matching well 
with profile 701 (discussed above). 

The laser power may then be changed, arid the measurements above 
repeated to measure a junction depth at a different active doping concentration. In 
this manner, the doping profile is determined as a plot of the results of the 

2 0 meastirements at different laser powers. 

In the second embodiment, operation 246 is performed by the acts 
illustrated in FIG. 6. Specifically, two interference signals are simultaneously 
me£isured. In act 621 the output of simi detector 522a (see FIG. 5) provides a 
measwe of the reference beam phase with respect to front surface 153 of the wafer 
25 (see FIG. IE). Moreover, in act 622 the output of the lock-in amplifier 525 (see 

FIG. 5) gives the difference signal based on the interference between the reflection 
1 63 by excess carriers at the junction and the reference beam 531. 

In act 623, both outputs are plotted on the same graph (e.g. FIG. 7B), and 
the phase shift between the two outputs is measured. In act 624, the phase shift is 

3 0 converted to a junction depth based on knowledge of the laser wavelength and the 

semiconductor index of refraction according to the formula 
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_ (phase _ shift) 



' 2kn 



where n is the index of refraction of the semiconductor and k=(2n)/X , where A, is 
the wavelength of probe beam 152. 

Ii» act 625, reference arm phase is then set to the position of maximum 
5 signal in detector 522a, corresponding to the reference beam path length being 
equal to the surface reflection path length. This makes the measurement output 
equivalent to that obtained using the front surface reflection method (first 
embodiment). In operation 626, one or more acts (described above in reference to 
the first embodiment) are used to verify the junction depth, and to determine the 
1 0 best fit to a doping profile. 

Line 701 (FIG. 7A) illustrates a typical doping profile, measured using SRP 
fi-om an ion implant under the conditions of 500 eV Boron ions at a dose of 1 x 
10'^ ions/cm^ annealed 10 seconds at 1000°C. Curves 702A-702C show the 
sunulated excess carrier concentration /cm^ for profile 701. with respective laser 

15 powers of 5 mW (for curve 702A), 20 mW (for curve 702B) and 50 mW (for curve 
702C), witii a spot diameter of 2 ^m. Applying eqn. (16) to the profiles for the 
excess carrier distributions at 5, 20 and 50 mW provides expected signals as a 
fiinction of reference arm phase of curves 712A, 712B. and 712C respectively. 
Also shown is curve 71 1, the cosine of the reference arm phase, corresponding to 

2 0 equation (1 9), with the zero point defined as Zref=Zs. 

A phase shift of 8=73° is observed in FIG. 7B. Solving using the above 
defined procedure. 



2 _ (phase _ shift) 
' 2kn 



where n=3.42 is the index of refi-action of silicon and k=(2n)/X , where the probe 
beam wavelength X = 1 .48 ^m, providing the above-mentioned value of 440 
angstroms. 
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The advantage of the second embodiment over the front-reflection method 
(first embodiment) is that it provides an absolute measure of the active carrier 
concentration depth corresponding to a set generation laser power level in terms of 
the phase shift of the cosine-shaped signal in the plot of signal versus reference arm 
5 length. Its disadvantages are the added complexity of the reference arm and the 
necessity of using a laser with a coherence length longer than the difference in 
physical length between the two arms of the interferometer, those being the 
reference arm and the measurement arm. 

Note that although in one specific embodiment, a reference beam having a 

1 0 variable phase has been described as being generated by an interferometer having 
two arms, any other device that can generate a variable phase reference beam can 
be used in accordance with the principles described herein. 

As noted above, in the second embodiment, interference with a reference 
beam 53 1 having a variable phase is used, and 50% of the light in probe beam 152 

15 is redirected towards competisator 506 and mirror 513 thereby to form the 

reference beam arm of the interferometer. Compensator 506 is used to set the 
polarization of reference beam 531 with respect to the remaining portion of probe 
beam 152 that is redirected toward objective lens 515 and wafer 516. Mirror 513 is 
mounted on a piezoelectric element 504, thereby to allow the length of the 

2 0 reference beam path to be adjusted electronically over the range of at least a 
wavelength (or several wavelengths if the signal is to be averaged over several 
cycles). 

The polarization of reference beam 531 is set orthogonal to probe beam 152 
(by compensator 506). Probe beam 1 52 and reference beam 53 1 interfere in 

2 5 polarizing beam splitter 52 1 , which is oriented at 45° v/ith respect to the 

polarization of the two beams 532 and 53 1 (wherein beam 532 represents the arm 
of the interferometer including the path to wafer front surface 153). This provides 
sum and difference beams at detectors 523a and 523b. The detector currents, 
which are proportional to the powers at the detectors, are converted to voltages in 

3 0 transimpedance amplifiers 523a and 523b. Amplifier 524 takes the difference in 

the voltages from amplifiers 523a and 523b. The difference signal is fed to lock-in 
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amplifier 525 thereby to generate the interference signal being supplied to 
processor 527. 

The physical principles that relate to the methods described above are as 
follows. The electric field at each of detectors 522a and 522b is 



E = E,+E,±E,, 



where Es is the portion 162 of probe beam 152 reflected from front surface 153 
(not including the excess carrier concentration at the surface), and Ej is the portion 

10 163 of probe beam 1 52 reflected by excess carriers m the junction (including the 
excess carriers near the surface), and Eref is the component from the variable phase 
reference beam, which is zero if there is no variable phase reference beam. 

The ± in equation 1 comes from the polarizing beam splitter 521 (FIG. 5) 
that provides sum and difference components to the two detectors 522a and 522b. 

15 The power in these detectors is 

p*={e3+e,±e«,Xe;+e;±e;^) 

Multiplying out equation 2, the signals m the sum and difference detectors 522a 
20 and 522b are 

P* =|E.r HEif -^lE-r +(e.e; ^e:e,)±(e.e;„ )±(e,e;. +e;e„) 

In the first embodiment, the surface reflection 162 is interfered with the reflection 
25 1 63 by the excess carriers, and hence the reference field is Eref is zero. Therefore, 
the squared amplitude of the reflection 163 from the junction is negligible, and the 
squared amplitude of the reflection 162 from fix)nt surface 153 does not appear at 
the modulation frequency. Thus, the only component seen in lock-in amplifier 525 
(FIG. 5) is 



(1) 



(2) 



(3) 
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P3.^=(e3e;4-e;e^) 



(4a) 



In the second embodiment, a reflection of the reference beam 53 1 by mirror 5 13 is 
interfered with the reflection 163 from the junction, the signal at the lock-in 
5 amplifier 525 is the difference between the signals from the plus and minus 
detector (e.g. from detectors 522a and 522b). Taking the difference, and 
recognizing that the i^^E^^^ + E[E^^^ ) term does not appear at the modulation 
frequency the resulting signal is 

10 P^,_j=P*-P-=2(EjE;,+E;E„,) (Sa) 

Finally, the junction term may be considered as a continuous distribution. 
Altemately, it may be considered to be composed of a surface and profile term, 
yielding 



PsH=[(Ejs+Eip)E;+(E;3^-E;pK] 



(4b) 



Pfef-i == 2[(Ej, + Ejp + (e;^ + e; ] (5b) 

2 0 Ignoring the time dependent part of propagation, the incident light from probe 
beam 152 has a phase at the surface of 



E- = E e^*^* 



2 5 where the amplitude is Eo= VP , with P the power of probe beam 1 52 at the wafer 
surface, and the wave number is k=27t/X, with X the wavelength of probe beam 
152. The reflected electric field from the boundary between air and the silicon 
surface 153 (FIG. IE) is 



(6) 
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where the amplitude is Eo= Vp , with P the power of probe beam 1 52 at wafer 
surface 153, nso is the index of refraction of silicon, Zs is the path length to the 
silicon surface, and the wave number is k=2n/X, with X the wavelength of the probe 
beam 152. 

As noted above, excess carrier profile 164 (FIG. IB) may be modeled as a 
set of thin layers of excess carriers, such as layers 164A-164T (FIG. IE). The 
reflection firom interface 165J between the j* layer 164J and (j+l)"' layer 164J+I is 



n +n 9n ^' (s) 



m 



where the index of refraction of a layer is nj =nso+An3 , where the change i 
index due to the excess carrier concentration is An, = ySN, , where is the 
15 excess carrier concentration in layer 1 64 J and the factor p is given by 



where q=I .602x10-'' is the electron charge, N is the carrier concentration per cm^ 
20 (the factor of 10^ converts to per m^ to allow use of MKS units), 80=8.86x10"" 
Farads/meter is the dielectric constant of free space, 83=11.7 is the relative 
dielectric constant of silicon, m*=5xl0-" kg is the effective carrier mass, and 
co=27cc/7, is the radial frequency, with c=3xl0'° cm/sec the speed of light. Note 
that concentration N is same as concentration nethat has been described above. 
2 5 The reflected electric field is the siun from all layers 1 64 A- 1 64T, 



(10) 
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where Zs is the distance to the surface and tg = 2^^^y^^ ^ ys the surface 

transmission. In the limit as the thickness of layer 164J approaches zero, applying 
8 into 10, 



(11) 



If the surface excess carrier concentration rises from zero to the.Ns over a distance 
Zf that is short compared to the wavelength, the integral in equation 1 1 can be 
10 broken into two parts. 



Ei=-(E.e'^-)^ti 



■so 



,j2knz 



|dZ+ J. 

2« 



dz 



(12) 



15 



which reduces to 



E,=-(Eoe'=^)J-« 



^ dz 



(13) 



20 



There are two terms in the above equation 13. One is from the near-surface, and 
depends only on the excess carrier concentration at the surface Nj. The second is 
the Fourier transform of the derivative of the excess carrier profile N. The electric 
field amplitude of the reference beam 531 is 



(14) 



25 where Zref is the length of the reference arm. Substitution into the above equations 
gives the signal power for interference between the surface and the excess carrier 
profile as 
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N3 + fcos(2knz)— dz 
o dz 



(15) 



10 



15 



20 



25 



and for interference between the reference arm and the excess carrier profile, 



"so 



N3 cos(2kz,,,)+/cos[2k(nz-z,,)]^dz 



Comparing equations (15) and (16), it is seen that when Zref=0, the two 
equations are identical. This condition is met when the path length for the 
reference beam is equal to the path length to the front surface 153 (see FIG. IE). 
This may be measured by simultaneously monitoring the output of either detector, 
523a or 523b. This output gives a signal that is the interference between the 
reference beam and the surface reflection. This term is the second term from the 
right in equation (3), and takes the form 



(E3E;, +E;E^)=2Pcos[2k(z, -z^)] 



(17) 



(16) 



When the reference and surface reflection paths are equal, the analysis to determine 
the semiconductor properties is identical for both embodiments. 

The above-described methods can be used for measuring non-dopant 
amorphizing implants (defined to be implants of ions that cause damage but do not 
dope silicon). Examples of non-dopant amorphizing implants are implants that 
cause damage when the high energy ions hit the crystal and stop. The ions have 
energies of several thousand eV, and the chemical bonds in the crystal are a few 
eV. Hence, the high energy ions are able to break the crystal lattice bonds, causing 
damage of silicon or germanium. 

The purpose of these implants is as follows. To form a shallow ion 
implanted layer with a dopant atom such as boron, it is essential to confine the 
boron ion implant to a very shallow surface layer (e.g. only a few angstroms thick). 
Typically, boron atoms are accelerated to a very low potential (e.g. a few hundred 
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eV) so that they stop within a few angstroms of the surface. However, because 
silicon has a crystal structure, certain directions present "channels" - long open 
paths through the crystal (typically the length of the crystal structure since channels 
are an inherent property of the crystal) in the space between the ordered silicon 
5 atoms. Some boron atoms will scatter off silicon atoms as they penetrate the 

silicon surface and follow the channels to considerably greater depths than desired. 

To prevent channeling, a non-dopant implant - often precedes the boron 
implant. The non-dopant implant creates enough damage to destroy the crystal 
structure, forming an amorphous layer that is characterized by the amorphization 

1 0 depth (a depth to which the crystal structure has been destroyed by ion implant 

damage and replaced by the amorphous layer). Amorphization depth is difficult to 
measure by existing means. The typical method for measuring amorphization 
depth is to create a very thin cross-sectional slice and look at it with transmission 
electron microscopy. This is a slow, tedious and destructive procedure 

15 The amorphous layer may be modeled as a high bandgap layer. As stated 

by Sze, "Physics of Semiconductor Devices," page 827, "The difference between 
crystalline and amorphous Si is dramatic; the former has an indirect bandgap of 1.1 
eV, whereas hydrogenated a-Si has an optical absorption characteristic that 
resembles the characteristic expected for a crystal with a direct bandgap of 1 .6 eV." 

2 0 Prior art FIG. 8 A (see figure 34 on page 828 of Sze), shows that at the preferred 

generation laser wavelength, 830 nm, the absorption in the amorphous silicon (line 
801) is only about 10% of that in the silicon (line 802). Thus, the amorphous 
silicon may be modeled as a transparent layer in one embodiment. 

Substituting into equation 15 for the signal in the case of interference with 

25 the surface (first embodiment). 
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where d is the thickness of the amorphous layer and the excess carrier 
concentration in the silicon under the amorphous layer is assumed to rise to a level 
AN in a distance very small compared to l/2knso, where n^o is the index of 
refraction of the silicon and k is the wave-number, equal to where "k is 

the wavelength. Here is it assumed that the index of refraction of the.silicon and 
the amorphized layer are approximately equal, so that n^o represents both layers. 

The above relation indicates that the signal will vary as the cosine of the 
thickness of the layer multiplied by 2kn. A calibration curve can therefore be 
established based on a measurement of reference samples to deteimine the signal as 
a function of depth. The curves in FIG 8B show examples of calibration curves for 
Ge (line 810) and Si Oine 81 1) [these designators have to be added to fig 8B; right 
now the lines are labeled "Si implant" and "Ge implant"] implants into siUcon. 
These curves are created by measuring the signal on reference samples created by 
ion implanting silicon or germanium into silicon. The amorphization depth can be 
calculated using the TRIM, program, available from J. F. Ziegler, Mail Stop 28- 
024, IBM-Research, Yorktown, NY, 10598. USA. For example, suppose a 
measurement is made on a sample that has been implanted with silicon. At a 
generation laser power of 90 mW the measured signal is 460 microvolts. This 
would indicate an amorphization depth of 320 angstroms. 

The above-described methods can be used to measure the active dopant 
profile in another special case, formed by the dose range of approximately 5x1 0'° 
to 5x1 ions/cm^ (often called the "low dose" range). Control of the dose of 
these implants is critical because they are used to adjust the turn-on voltage 
(threshold voltage) of field effect transistors. Small variations in these "low dose" 
implants can resuh in tum-on voltages that are out of the operatmg range of the 
integrated circuit. 

Applicants recognize that such "low dose" implants cannot be measured by 
thermal wave methods of the type described by Opsal (US Patent 4,854,710) with 
high sensitivity, because such methods rely on the decay of the propagation of 
thermal waves due to the damage resulting from the ion implant. However, low 
dose implants create relatively little damage, and Opsal's thermal wave methods 
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generally have very weak sensitivity in the low dose range (the sensitivity, defined 
as the percent change in signal divided by the percent change in dose, is typically 
about 0.2 for thermal wave methods; a sensitivity of 0.5 is usually considered the 
minimum usable level, and sensitivities > 1.0 are considered desirable). 

A significant reason for the weak sensitivity of OpsaPs thermal wave 
methods is that they are sensitive to mobility but not lifetime. This may be shown 
as follows. The carrier distribution is a solution to the time dependent diffusion 
equation, 

d^n n _ 1 an 

az' Dt " D at 



where n is the excess carrier concentration, D is the diffusion constant, 

withD = [^*^^yq^ ' where kb is Boltzmann's constant, T is the temperature, q the 

electron charge, and the mobility, t is the lifetime. For a periodically excited 
15 carrier concentration at a radial frequency co, the carrier concentration is 

n(2, t) = n(z)e^^ , where t is the time. This gives a diffusion equation of the form 



a'n ri .G)l ^ 



20 At high fi-equencies, (0»l/t (for an excess carrier lifetime of 10"^ seconds, this 
represents a frequency greater than about 15 kHz). The solution is of the form 
n(z, t) = n^e^^*"*"*^^ , where k^=a)/D . This is a wave propagating solution 

whose propagation constant is a function of the mobility, since D = ^'^^'^q^ • 

This is the region where the thermal wave method operates, since that method 
25 relies on a propagating wave. 
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Conversely, at low frequencies, the solution is of the form 
n(z. t) = Hoe ^/D7gj(«.i) y^jg gj^Q^g ^ gj^^j^ gp^^-^j variation that is a fiinction of 
both mobility and lifetime. The import of this result is that the lifetime of a 
semiconductor is several orders of magnitude more sensitive to defect density than 
the mobility, so a measurement sensitive to lifetime will show greater sensitivity to 
the damage caused by low doses of ion implantation than a measurement sensitive 
to the mobility. 

Therefore, there is a transition from a wave behavior to a diffusion behavior 
(that is sensitive to carrier lifetime) when co = 1/x , i.e. when f = (I/27it). 
Therefore, in one embodiment, the modulation frequency is preselected to be any 
frequency in conformance with the formvila f<. (1/27it), 

A graph illustrated in FIG. 9A shows the results of calculations based on a 
solution to the diffusion equation assuming the ion implanted layer at the surface 
has a reduced lifetime compared to the bulk silicon under the implanted region. 
15 The horizontal axis is the depth in microns and the vertical axis is the carrier 
concentration per cubic cm. Note that for small changes in lifetime in the 
implanted layer, the carrier concentration is nearly constant. However, as the 
lifetime in the implanted layer becomes shorter, the surface concentration starts to 
drop very quickly. Over this range, the surface concentration drops quickly. At a 
2 0 certain point - about lO"' ' seconds in this model - the carrier concentration 

becomes relatively independent of lifetime again. However, the point of reflection 
shifts from the surface to the boundary between the unplanted layer and the bulk 
silicon, hi this regime, the signal levels off. 

Another graph illustrated in FIG. 9B shows the signal for various low dose 
2 5 implants, including a variety of species (B - Boron, As - Arsenic, P - 

Phosphorous, and BF2 - Boron Fluoride) as a fiinction of dose. For the lower 
doses, the signal is seen to drop relatively rapidly. For the higher doses, the signal 
flattens out This is consistent with the above model. The signal in the steep part 
comes from the surface reflection - the first term in the parenthesis in equation 1 5. 
30 As the signal flattens and begins to come from reflection from the interface 
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between the implant and the bulk, the interference term contributes. This explains 
why the signal changes phase, and is negative at the higher dose values. 

The graph in FIG. 9B may be used to calibrate the measurement in a 
manner similar to that for the amorphous measurement. The signal as a function of 
5 dose is measured and stored as a graph similar to FIG. 9B. For example, a signal 
of 1 00 for a Boron implant corresponds to a dose of 5e 1 L 

One implementation, software uses a power curve (see curve 404 in FIG. 
4C) obtained by measurement of an interference signal (amplitude and phase, 
wherein phase is used to determine a sign to be used with the amplitude) as a 

1 0 function of the power of generation beam 151. Specifically, the software compares 
the measured power curve to a number of preexisting power curves that have been 
obtained by simulation (see curves 404-409 in FIG. 4E) to find a match. The 
matched power curve obtained by simulation determines the junction depth. The 
preexisting power curves are obtained as described above in reference to FIGs. 4D 

15 and4E. 

Numerous modifications and adaptations of the above-described 
embodiments will become apparent to a person skilled in the art of semiconductor 
physics. For example, although computer 103C is described as being programmed 
with one or more specific equations, computer 103C can be programmed with 

2 0 other equations described herein, or with one or more equations that approximate 
any of the relations between material properties as described herein, for use with 
measurements performed by profiler 103 while creating a diffusive modulation of 
charge carriers in a wafer under measurement. For example, an approximate 
equation used by profiler 103 to measure a material property can be obtained by 

2 5 curve-fitting to measurement data from reference wafers, or by curve-fitting to data 
obtained from a numerical model, or both depending on the specific 
implementation. 

Therefore, numerous such modifications and adaptations of the above- 
described embodiments are encompassed by the attached claims. 
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Claims 

What is claimed is 

1 . A method for performing a measurement in a region of a wafer 
having a plurality of background carriers, the method comprising: 

creating a plurality of excess carriers in the region, the number of 
excess carriers in the plurality being modulated in time at a predetermined 
frequency that is sufficiently small to ensure that a majority of carriers 
move out of the region by diffusion; and 

measuring amplitude and phase of an interference signal, the 
interference signal being obtained by interference between: 
a reference beam; and 

a portion of a probe beam of electromagnetic radiation 
reflected by the plurality of excess carriers in the region, the portion 
of the probe beam being modulated in phase with modulation of the 
plurality of excess carriers. 

2. The method of Claim 1 wherein: 
the reference beam is formed by another portion of the probe beam 

that is reflected from a front surface of said substrate. 

3. The method of Claim 2 fiirther comprising: 
usmg predetermined data to look up a junction depth corresponding 

to measured values of the amplitude and phase obtained by the "measuring 
amplitude and phase"; 

wherein the predetermined data relates the measured value to a 
known junction depth of a predetermined substrate. 

4. The method of Claim 3 further comprising: 

generating the predetermined data by operating a simulator to 

determine a profile of excess carriers as a function of depth for a given 
dopant profile in the substrate. 
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5. The method of Claim 4 wherein said generating further comprises: 
multiplying: 

a derivative of the profile of excess carriers with 
5 respect to depth from the front surface; and 

cos (2knz), wherein z is the depth, k=27t/X, X is the 
wavelength of the probe beam, and n is the index of 
refraction of the substrate; 
integrating the product of the multiplying with respect to depth from 
10 the front surface; and 

multiplying a result of the integrating with a constant to determine a 
simulated value of the amplitude and phase of the interference signal. 

6. The method of Claun 1 fiirther comprising: 

15 repeating the measuring at a plurality of times during the 

modulation; and 

computing an average value of a property of the region from a 
plurality of measured values of the property obtained during the repeating, 
and during the measuring. 

20 

7. The method of Claim 6 fiirther comprising: 

using predetermined data to look up a junction depth corresponding 
to the measured value. 

25 8. The method of Claim 1 wherein the region is henceforth referred to 

as first region, the method fiirther comprising: 

creating a plurality of excess carriers in a second region adjacent to 

the region; 

repeating the "measuring amplitude and phase" in the second 
30 region; and 
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computing a difference between a second measured value in the 
second region obtained during the repeating and a first measured value in 
the first region obtained during the "measuring amplitude and phase." 

5 9. The method of Claim 8 further comprising: 

comparing the difference with a predetermined limit; and 
changing a process parameter used in fabricating the substrate if the 
difference exceeds a predetermined limit. 

10 10- The method of Claim 1 further comprising: 

annealing the substrate to activate dopants; 
wherein the measuring is performed after the annealing. 

11. The method of Claim 1 further comprising: 

1^ changing an average number of excess carriers in the plurality; 

repeating the "measuring amplitude and phase"; and 
using a plurality of predetermined curves fitted to the predetermined 
data to look up a junction depth that best fits a first measured value 
determined by the "measuring amplitude and phase" and a second measured 
20 value determined by the repeating; 

wherein each predetermined curve relates the amplitude and phase 
to one of a plurality of known junction depths of a predetermined substrate. 

12. The method of Claim 1 wherein the predetermined frequency is in 

2 5 conformance with the formula: 

f± (1/27CT) 

where f is the fi-equency, and x is the lifetime of a carrier in the substrate. 

13. The method of Claim 1 wherein said plurality of excess carriers are 

3 0 created by a generation beam formed by a laser, the generation beam having a first 
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wavelength A.g and the probe beam havinjg a second wavelength Xp, the second 
wavelength Xp being in conformance with the formula: 

> [ (1 0 apPp;ip)/( agPg) ] [wg/wp]' 
wherein ap and ag are the absorption coefficients in the substrate of the 
5 probe beam and the generation beam respectively, Pp and Pg are the powers of the 
probe beam and the generation beam respectively, and Wg and Wp are radii of focal 
spots of the generation beam and the probe beam at a front surface of said wafer. 

14. The method of Claim 13 wherein: 

10 the generation beam includes a plurality of photons that have energy greater 

than the bandgap energy of the wafer. 

15. The method of Claim 13 wherein: 

the generation beam includes a plurality of photons that have energy greater 
15 than the bandgap energy of crystalline silicon but less than the bandgap energy of 
amorphous silicin. 

16. The method of Claim 1 wherein: 

the probe beam beam includes a plurality of photons that have energy lower 
2 0 than the bandgap energy of the wafer. 

17. The method of Claim 1 wherein: 

the probe beam beam includes a plurality of photons that have energy 
approximately equal to the bandgap energy of the wafer. 

25 

1 8. The method of Claim 1 wherein: 
wherein said frequency is in accordance with the formula: 

f^ (l/27tT) 

where f is the frequency, and x is the lifetime of one of the excess carriers. 

30 
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1 9. The method of Claim 1 wherein: 
the probe beam is coherent; 

the reference beam is coherent with respect to the probe beam; and 
the reference beam has a phase that is variable independent of the phase of 
the probe beam. 

20. The method of Claim 1 9 further comprising: 
changing the phase of the reference beam; and 
repeating the measuring. 

2 1 . The method of Claim 20 wherein: 

the measuring includes using a phase detector to detect the difference in 
phase between: 

a first interference signal obtained by interference of the reference 
beam and another portion of the probe beam reflected by the front surface; 
and 

a second interference signal obtained by interference of the 
reference beam and the portion of the probe beam reflected by the plurality 
of excess carriers. 

22. A method for perfonning a measurement in a region of a wafer 
having a plurality of charge carriers, the plurality consisting of a first number of 
charge carriers, the method comprising: 

stimulating the region, the plurality consisting of a second number 
of charge carriers during the stimulating; 

using a measurement device to obtain a measured value of a signal 
affected by the stimulating; 

operating a simulator to generate a simulated value of the signal 
based on: 

(i) a predetermined profile of concentration of active 
dopants in the region; and 
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(ii) conditions of the stimulating; and 
comparing the simulated value with the measured value to 
determine if the simulated value matches the measured value. 

5 23. The method of Claim 22 further comprismg: 

repeating the operating for a set of predetermined profiles including 
the predetermined profile to jgenerate a set of simulated values including the 
simulated value; 

repeating the comparing with each simulated value in the set to 
1 0 determine the simulated value closest to the measured value to be a match. 

24. The method of Claim 22 further comprising: 

repeating the operating with a number of additional predetermined 
profiles until an additional simulated value and the measured value differ 
15 by less than a predetermined amoimt; and 

determining the additional simulated value to be a match. 

25. The method of Claim 22 further comprising: 
computing a value of a property of the region based on the 

2 0 predetermined profile. 

26. The method of Claim 22 further comprising: 

driving a control signal to a device selected from the group 
consisting of a wafer processing imit and an annealer based on a result of 
2 5 said comparing. 

27. The method of Claim 22 wherein: 

the measurement device includes an interferometer; and 
the signal is an interference signal obtained by interference between: 
30 a reference beam; and 
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a portion of a probe beam of electromagnetic radiation 
reflected by a plurality of excess carriers in the region, the plurality 
of excess carriers being created by the stimulating. 

28. The method of Claim 27 wherein: 

the reference beam is formed by another portion of the probe beam 
that is reflected from a front surface of said substrate. 

29. The method of Claim 22 wherein: 

the stimulating includes creating a plurality of excess carriers in the 
region, the number of excess carriers in the plurality being modulated in 
time at a predetermined frequency that is sufficientiy small to ensure that a 
majority of carriers move out of the region by diffusion. 

30. The method of Claim 29 wherein: 

the operating includes determining a profile of excess carriers as a 
function of depth in the substrate. 



3 1 . The method of Claim 30 further comprising: 

using the profile of excess carriers to determine tiie simulated value. 

32. An apparatus for performing a measurement in a region of a wafer 
having a plurality of background carriers, said apparatus comprising: 

means for creating a plurality of excess carriers in a region of the 
substrate, tiie number of excess carriers in the plurality being modulated at 
a frequency tiiat is sufficientiy small to cause a majority of carriers moving 
out of the region to transfer by diffusion; 

a source of a probe beam of elecfromagnetic radiation; and 
an interferometer located in a path of a signal obtained by 
interference between a reference beam and a portion of the probe beam 
reflected by tiie plurality of excess carriers in tiie region, tiie portion of tiie 
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probe beam being modulated in phase with modulation of the plurality of 
excess carriers. 

33. The apparatus of Claim 32 wherein the interferometer includes a 
lock-in amplifier that detects amplitude and phase of a signal obtained by 
interference between: 

the portion of the probe beam reflected by the pluralit>^ of excess carriers; 

and 

another portion of the probe beam reflected by a front surface of the 
substrate. 

34. The apparatus of Claim 32 further comprising: 
a detector coupled to the interferometer, wherein the detector is used to 

measure a difference in phase between: 

a first interference signal obtained by interference of the reference 
beam and another portion of probe beam reflected by a firont surface of the 
substrate; and 

a second interference signal obtained by interference of the 
reference beam and the portion of the probe beam reflected by the plurality 
of excess carriers. 

35. The apparatus of Claim 32 fiirther comprising: 
a computer coupled to the interferometer, the computer being programmed 

to use predetermined data to look up a junction depth that corresponds to a power 
25 of the portion of probe beam measured by the interferometer. 

36. The apparatus of Claim 32 wherein said frequency is in accordance 
with the formula: 

f <. (10/27CT) 

3 0 where f is the frequency, and x is the lifetime of one of the exccess carriers. 



5 



10 



15 



20 



-59- 



wo 00/57159 




PCT/USOO/07357 



37. The apparatus of Claim 32 wherein: 

said means for creating includes a laser that generates a generation beam of 
electromagnetic radiation having a first wavelength Xg; and 

the probe beam has a second wavelength Xp, the second wavelength Xp 
5 being in accordance with the formula: 

> [ (10 apPp>.p)/( agPg) ] [wg/wp]2 
wherein ap and ag are the absorption coefficients in the substrate of the 
probe beam and the generation beam respectively, Pp and Pg are the powers of the 
probe beam and the generation beam respectively, and Wg and Wp are radii of focal 
0 spots of the generation beam and the probe beam at a front surface of said wafer. 



CXJO: <WO 00571S9A1_I.> 



-60- 



wo 00/57159 ^pp'/USOO/07357 

1/14 



101 



104 



r 

100 



WAFER PROCESSING 



UNIT I 

1011^ ION 
IIMPLANTER 



107— 



105 



105/106 



RAPID 
THERMAL 
ANNEALER 



10^ ^ 
106 



—108 



r 
I 



103 
4. 



ACTIVE DOPANT PROFILER 



PROGRAMMEDI 

COMPUTER ^-103C 
_i 



109 



CONTROL SIGNALS 




103M 



FIG. lA 



suBSTmrre sheet (rule 26) 



wo 00/57159 _ PCT/USOO/07357 



2/14 



r 



operate a simulator to generate a simulated value 
of a to-be-measured signal based on a predetermined 

profile of the concentration of active dopants in a 
region of a semiconductor substrate, and based on the 
conditions to be present during stimulation of the region 

111 



113 



change predetermined profile 
used in operating the simulator 



110 



that 



USE SIMULATION 
TO GENERATE 
^REDETERMINED DATA 



done with a set of" 
predetermined profiles 
the semiconductor substrate 
is likely to have)? 



125 



stimulate a region of semiconductor 
substrate to change the number of 
charge corners in the region 



r 



Prepare 
Semi- 
conductor 
Substrate 



-126 



use a measurement device to obtain a 
measured value of a signal that is affected 
by the stimulotinq 



-141 



compare the measured value with one 
or more of the simulated values to identify 
the closest simulated value 



determine a value of a property of the 
region based on the predetermined profile 
that generated the closest simulated value 



142 



'143 



MEASURE A 
PROPERTY THAT 
IS AFFECTED BY 

STIMULATION 

140 




drive a control signal to a wafer 
processing unit or an annealer 



■161 



FIG. IB 



)OCIO: <WO_0057159A1„I_> 



SUBSTITUTE SHEET (RULE 26) 



wo 00/57159 



/USOO/07357 



3/14 



prepare semiconductor substrote (by adding 
dopants in wafer processing unit and optionally 
by heot treatment in on onnealer) 



'125 



stimulate a region of semiconductor substrate - 
to change the number of charge carriers in the 
region (from a first number to o second number) 



^126 



use 0 measurement device to obtain a 
measured value of a signal that is affected by 
the stimulating (e.g. indicative of the difference 

between the first and second numbers) 



'171 



USE SIMULATOR 
TO MEASURE A 
PROPERTY AFFECTED 
BY STIMULATION 

170 



operate a simulator to generate a simulated 
value of the signal based on a predetermined 
profile of the concentration of active dopants in 
the region, and based on the conditions present 
during the stimulating 



'172 



^•173 

does 
^simulated value^ 
match measured 
value? 



174 



change predetermined profile 
used in operating the simulator 



compute a value of a property of the region 
based on the predetermined profile that 
produced the matching simulated value 



'175 



^160 



Yes 



does property 
value match 
specifications?^ 



drive a control signal to a wafer 
processing unit or oh onnealer 



161 



FIG. IC 



SUBSTITUTE SHEET (RULE 26) 



wo 00/57159 



4/14 



PCT/USOO/07357 



152 



153 
130 



155ZZ.151 

y^l05/106 



156 



120 



FIG. IE 



152- 



z"iF:E:"i35r:EEE?i;;;;;:EEE=E 



162 



153 



156 



FIG. IF 





0 

zj- 

2 



FIG. ID 



153 



-164 



130 



156 



FIG. IG 



OCID: <WO_0057159A1_L> 



SUBSTITUTE SHEET (RULE 25> 



wo 00/57159 



5/14 



'AJSOO/07357 



201 Start 



Use simulation to --210 
generate predetermined 
dato based on profiles 



r 



205 



221 ' Fornn doped regions 



222 ' Anneal the wafer 



L. 



223 



Prepare 
Semiconductor 
Substrate 



± 



224-n- Position wafer relative to 
active dopant profiler 



3: 



220 



J 



230- 



Create in a region of the 
wafer, carriers in a concentration 
that is modulated at a frequency 
selected to be low enough to 
avoid the creation of a wave 



240 



242 



244 
J— 



Generate & Focus a probe beam 
— - (e.g. photons having energy 
lower than bandgap 
energy of semiconductor 
material in the region) 



Change a parameter 
used in generating 

reference beam 
(e.g. change phase) 



T 



1 



I 



Measure Property 
Affected By The 
Created Charge Carriers 



241 



Measure phase and intensity 
of interference between reflection 
of probe beam by excess carriers, 
and another portion of probe 
beam (colled reference beam) 



I 



245 



Change a parameter 
used in creating 

(e.g. overage carrier 
concentration or 

locat ion of region) 

t 



243 



Use one or more 
measurements to look up 
semiconductor property 
(e.g. junction depth/profile) 
from predetermined doto 



L_! - — — ^ 



-246 



200 



-260 

/ Do ^' ^ 
'measurements(s) 
Indicate tliat wafer is withir 
specifications^ 



FIG. 2A 



Reject wafer and 
adjust adding of dopants 
or annealing or both if 
necessary 



-261 



Identify wafer as 
acceptable and 
process wafer 
further in the 
normal manner 

C — 

262 



263' 



SUBSTITUTE SHEET (RULE 26) 



6/14 



PCT/USOO/07357 



Generate a beam of photons having 
energy greater than the bandgap 
energy of semiconductor moteriol 



t 
231 



—232 



Modulate 
bandgap b( 
selected 


the above 
3am at the 
frequency 






Focus the above bandgap' 
beam on the region 



■233 



•230 



FIG. 2B 



Transform the 

interference 
measurement(s) 



i 



-251 



252 



find the location of transformed 
measurement on a curve representing 
a correlation between transformed 
power and junction depth for a 
selected concentration of excess carriers 



Read the junction 
depth from location 



■253 



■246 



FIG. 2C 



SUBSTITUTE SHEET (RULE 26) 



wo 00/57159 



7US00/073S7 



7/14 



100 

SA 

50 
SI 

Siqnol 0, 

-50 

-100 
SN 

-150 



Signal vs. depth @10 mW 



- 255A^^^X^^^^ y=-3016x + 130.88 


^'^"^'^^ 


^2551 R2 = 0.9957 


.00 0.01 0.02 0.03 0. 


04'*^S^:^'06 0.07 0. 











Depth (um) 



-255N 



FIG. 2D 



500 
400 
300 

Signal o 
(uV) -100 
-200 
-300- 
-400 
-500 



Signal vs. depth @ 40 mW 









v=-1442bx + /32.4 






-100 0.01 0.02 0.03 0. 






Q'l 0.q8^»^^6 q!07 0. 















Depth (um) ^2 



FIG. 2E 



S3- 



1500- 

1000- 



25y^ Signal vs. depth @ 100 nnW 



Signal 
(uV) 



500' 



0 
C 

-500- 
-1000- 



^25T~~"^*^^^^ Slope = 35,400 uV/um 




y=-35378x + 2026.7 








).00 0.01 0.02 0.03 0. 


04 0.05 o^ot^O^ 0 




^V3 


Depth 


(um) 



08 



•257N 

FIG. 2F 



SUBSTITUTE SHEET (RULE 26) 



wo 00/57159 



8/14 



# 



PCT/USOO/07357 



310 



Receive a dopant profile 
(from simulation or from 
measurement (e.g. SIMS or SRP)] 
for a selected junction depth 



Use a simulator to determine 
a profile of excess carriers 
for the dopant profile and 

for selected concentration of 
excess carriers 



^311 



312 



Multiply (a) derivative of profile 
of excess carriers with respect to 

depth from front surface and 
(b) cos(2kn2) wherein z is depth, 
K= 2u/x,Xis wavelength of probe 

beam and n is refraction index 



313 



integrate product of multiplying 
with respect to depth from front surface 
and add excess concentration at 
surface to integral 



314 



Multiply (a) result of integrating 
and (b) constant based on 
calibration of measurement 
system to obtain an expected 
value of interference measurement 



'315 




FIG. 3 



SUBSTITUTE SHEET (RULE 26) 

OCIO: <WO_0057159A1_I_> 



wo 00/57159 



Log 
cone 



21 
20 
19 
18 
17 
16 



15 
0.00 



9/14 



Doping Profile 




0.02 0.04 0.06 
Depth (urn) 



AJSOO/07357 



0.08 



0.10 



FIG. 4A 




Signal 



(w) 



2.0E-04 
1.5E-04 
1.0E-04 



5.0E-05 
O.OE+00 
-5.0E-05 



404J 



404P 



Signol vs. power 



404A- 




^ 10 20 30 40 50 60 70 80 90 ibo 



Power @ wafer (mW) 



FIG. 40 



SUBSTITUTE SHEET (RULE 26) 



wo 00/57159 



10/14 



PCT/USOO/07357 





21- 



Doping Profile 



Log 
cone 
(/cm'^3) 



20 

19- 
18- 

17- 
IB- 



IS- 
0.00 




0.02 0.04 0.06 
Depth (um) 



FIG. 4F 



0.08 



0.10 



.0057159A1J_> 



SUBSTITUTE SHEET (RULE 26) 




Measure sum detector 
signal to determine 

phase and intensity of 
interference between 

front surface reflection 
ond reference beam 



Measure difference 

signal to determine 
phase and intensity of 

interference between 
excess carrier reflection 

and reference beam 



Determine phase shift between front surface- 
reference beam interference and reflection from 
junction-reference beam interference 



Convert the phase shift to a junction depth 



Set reference arm to zero phase position 



Conduct measurements as described for front - 
surface embodiment to verify junction depth and 
determine best fit to a doping profile 



'623 



■624 



625 



626 



'246 



FIG. 6 



SUBSTITUTE SHEET (RULE 26) 



wo 00/57159 



12/14 



PCT/US0O/O73S7 



< 

E 

o 



o 
o 



o 



excess holes -r electrons 



702C 




50 



^702B 9^ 



0.04 0.06 
Depth (urn) 

FIG. 7A 



0.08 



21 
20 
19 
18 
17 



-16 



15 



0.10 



4.00E-04 
3.00E-04 
2.00E-04 
^ 1.00E-04 
o O.OOE+00 
^ -1.00E-04 
-2.00E-04 
-3.00E-04 
--4.00E-04 



Signal vs. ref. phase 




Reference phase (deg) 



FIG. 7B 



X3CID: <WO 00571 59AT_J_> 



SUBSTITUTE SHEET (RULE 26) 



wo 00/57159 ^^;r/US00/07357 

13/14 



a-Si and Si absorption 



1.E+05 

T 

I l.E+04 

o 

8- l.E+03 
1.E+02 



O 
in 



0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 
Wavelength (um) 



•a-Si 
Si 



FIG. 8A 

(PRIOR ART) 



Depth (A) 
500 




Si implant^/^811 j^r^piant 

810 



400 500 600 700 800 
Signal (uV) @290 mW laser power 

FIG. 8B 



900 1000 



SUBSnrUTE sheet (rule 26) 



wo 00/57159 



14/14 



Cone. 
(/cm3) 




0 



0.1 



0.2 0.3 
Depth (urn) 



0.4 



0.5 



FIG. 9A 



signal 



10( 



-150 



r ■ . 

) — 


■ 




^ 












Mi 




•h 




1 





-*^B 9keV 
-*-B 20keV 
-^As 9keV 
P 20keV 
-»-As 20keV 
-^BF2 40keV 

FIG. 9C 



Nominal dose 

FIG. 9B 



SUBSTITUTE SHEET (RULE 26) 

OOCIO: <VVO_00S71S9A1_L> 



INTERNATION. 



CH REPORT 



In itlonai 



«ion No 
'357 



A. CLASSIFICATION OF SUBJECT MATTER 

IPC 7 G01N21/17. 



According to imemational Patent Classification f IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

IPC 7 GOIN 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched 



Electronic data base consulted during the international search (name of data base and. where practical, search terms used) 

WPI Data, PAJ, EPO-Internal . INSPEC 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Categoiy * Citation of document, with indicalian, where appropriate, oi the relevant passages 



Relevant to daim No. 



US 5 042 952 A (OPSAL JON ET AL) 
27 August 1991 (1991-08-27) 
column 5, line 7 - line 37 
claim 1; figures 1,2 

US 5 377 006 A (NAKATA TOSHIHIKO) 
27 December 1994 (1994-12-27) 
column 2, line 64 -column 5, line 2 

US 5 706 094 A (MARIS HUMPHREY J) 
6 January 1998 (1998-01-06) 
claims 1,5,8,35,38,41 



1,22 



1.22 



1,22 



□ 



Further documents are listed in the continuation of box C. 



ID 



Patent family members are listed in annex. 



° Special categories of cited documents : 

'A* document defining the general slate of the art which is not 
considered to be of particular relevance 

*E* earlier document but publistned on or after the international 
filing date 

"L" document which may throw doubts on priority c<atm(s) or 
which is cited to estatriish the publication date of another 
citation or other special reason (as specified) 

'O" document referring to an oral disclosure, use. exhibition or 
other means 

'P' document pubiished prior to the international filing date but 
later than the priority date ciaimed- 



T* later document pubiished after the international filirtg date 
or priority date and not in conflict with the application but 
cited to understar • ne principle or theory undertyir>g the 
invention 

"X* document of particular relevance; the claimed invention 
carmot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

"Y" document of particular relevance: the claimed invention 

cannot t>e considered to involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination beir^ obvious to a person skilled 
tn the art. 

document member of the same patent family 



Date of the actual completion of the intematior%al search 



4 July 2000 



Date of mailing of the international search report 



12/07/2000 



Name and mailing address of the ISA 

European Patent Office. P.B. 5816 Patentiaan 2 
NL - 2280 HV Rijswiik 
Tel. (431-70) 340-2040, Tx. 31 651 eponl. 
Fax: (+31-70) 340-O016 



Authorized officer 



Krametz, E 



Form PCT/ISA/210 (second shMU CJuly 1992) 



OCID: <WO_0057159A1J.> 



INTERNATl 



TIONAJL SEARCH REPORT 

lation on patent family membera 



Patent document 
cited in search report 



Publication 
date 



^tionat Application No 

us 00/07357 



Patent family 
inernDer(s) 



Put)lication 
date 



US 5042952 



US 5377006 



US 5706094 



27-08-1991 



27-12-1994 



06-01-1998 



US 


4636088 


A 


13-01- 


■1987 


DE 


3578824 


0 


30-08- 


-1990 


EP 


0162681 


A 


27-11- 


•1985 


TP 

Ur 


looZoZo 


C 


31-07- 


■1992 


jp 


3047704 


B 


22-07- 


■1991 


jp 


60256031 


A 


17-12- 


•1985 


EP 


0200301 


A 


05-11- 


•1985 


OP 


2002847 


C 


20-12- 


•1995 


JP 


61223543 


A 


04-10- 


•1986 


US 


4854710 


A 


08-08- 


•1989 


us 


4952063 


A 


28-08- 




OP 


4343045 


A 


30-11- 


•1992 


Jp 


3000729 


8 


17-01- 


2000 


JP 


5001988 


A 


08-01- 


■1993 


us 


5479259 A 


26-12- 


1995 


OP 


11511240 


T 


28-09- 


1999 


wo 


9708536 


A 


06-03- 


1997 


us 


6008906 


A 


28-12- 


1999 



Form PCT,1SA/210 (patent (am«y annex) (July 1992) 
0010: <WO 0057159A1_I_> 



